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We investigate the nature of the mid-IR excess for 31 intermediate-mass stars 
t:j- ■ that exhibit an 8 yum excess in either the Galactic Legacy Infrared Mid-Plane Sur- 

vey Extraordinaire or the Mid-Course Space Experiment using high resolution 
^ ■ optical spectra to identify stars surrounded by warm circumstellar dust. From 

CL|' these data we determine projected stellar rotational velocities and estimate stel- 

O ' lar effective temperatures for the sample. We estimate stellar ages from these 

^ ■ temperatures, parallactic distances, and evolutionary models. Using MIPS [24] 

, c^, ■ measurements and stellar parameters we determine the nature of the infrared ex- 

cess for 19 GLIMPSE stars. We find that 15 stars exhibit Ha emission and four 
^ ■ exhibit Ha absorption. Assuming that the mid-IR excesses arise in circumstel- 

^ ■ lar disks, we use the Ha fluxes to model and estimate the relative contributions 

00 . 

Q^ . of dust and free-free emission. Six stars exhibit Ha fluxes that imply free-free 
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emission can plausibly explain the infrared excess at [24]. These stars are can- 
didate classical Be stars. Nine stars exhibit Ha emission, but their Ha fluxes 
OO , are insufficient to explain the infrared excesses at [24], suggesting the presence 

of a circumstellar dust component. After the removal of the free- free component 



X 



in these sources, we determine probable disk dust temperatures of Idisk —300- 
800 K and fractional infrared luminosities of LinjL^ ~10^'^. These nine stars 
may be pre-main-sequence stars with transitional disks undergoing disk clear- 
ing. Three of the four sources showing Ha absorption exhibit circumstellar disk 
temperatures ~300-400 K, Lm/L^, ~10~^, IR colors -ft'-[24]< 3.3, and are warm 
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debris disk candidates. One of the four Ha absorption sources has i^-[24]> 3.3 
implying an optically thick outer disk and is a transition disk candidate. 

Subject headings: Line: profiles — Stars: Fundamental Parameters, Emission- 
line, Be, Rotation, Circumstellar matter — Techniques: spectroscopic 



Introduction 



With its launch in 1983, the Infrared Astronomical Satellite (IRAS) opened up a new 
wavelength window and enabled identification of many classes of stellar sources with circum- 
stellar material, providing new insights into star forma t ion an d evolution. In one of the most 
exciting discoveries enabled by IRAS, lAumann et al.l (Il984j ) reported the first detection of 
an infrared excess owing to dust in the Vega system, and this excess was later explained 
as a debris disk. To date, many debris disks have been found with IRAS, the Infrared 
Spac e Observatory (ISO), and the Syitzer Sp ace Telescope (Spitzer) (jBackman &: Paresce 
19931 : ISpangler et al.l 1200 ll : iMeyer et al.l |2004J ) . Most known debris disks show evidence of 
cool dust at wavelengths > 20 / im, and the majority of these surround intermediate-mass 
stars JBackman fc ParescelEggsl : Ichen et allbood : ImUenbrand et al.lboOSJ ). Only two IRAS 
Point Source Catalog debris disk sources, the intermediate-mass A-type stars (3 Pictoris and 



C Lep oris, exhibited a mid-IR (< 12 /xm) excess detectable by IRAS ( lAumann fc Probst 



199ll ). Spitzer observations of /i?A5'-discovered de bris disks indicate that most are devoid 
of warm dust and small particles (sizes < 10 /xm) (jChen et al.ll2006l ). The rarity of mid- 
IR excess sources may imply that dust particles emitting at terrestrial temperatures are 
uncommon around main-sequence stars. 

Advances in infrared technologies coupled with large-aperture ground-based telescopes 
and Spitzer' s mid-IR sensitivity with IRAC enable new searches for warm/hot circumstellar 
dust that IRAS did not detect. Even with such improvements in sensitivity over IRAS, 
only five known intermediate-mass stars (HD 3003, HD 38678, HD 75416, HD 105234, and 
HD 181296) are surrounded by debris disks with tem peratures analogous to our asteroid 



belt (> 200 K) (IChen fc Jural l200ll : iRhee et al.ll2007l and references therein). Spitzer ob 



servations have also revealed a number of mid-IR excess objects that a re primordial disks 



undergoing disk cl earing processes similar to the solar-mass star TW Hya (ICalvet et al.ll2005 



Furlan et al.ll2007l ). "Transition disks" may be evolved protoplanetary di sks with inner hole s 
caused by clearing of the initial gas-dominated disks from the inside out (IStrom et al.lll989l ). 
These disks would be optically thin at near- and mid-IR wavelengths a nd optically thick in 



the far -IR. These sources would show weak to no accretion signatures (jSicilia-Aguilar et al. 
2006bl ). Intermediate- mass candidate transition disks have only recently been discovered 
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( iHernandez et al.ll2006l : ISicilia-Aguilar et al.ll2006al : iHernandez et al.ll2007l ). Such transition 
disks are especially interesting as possible evolutionary links between the optically thick 
Herbig AeBe disk systems and the older optically thin debris disk systems. 



In lUzpen et al.l (l2005l . 120071 ) we exploited t his new mid-IR cap abihty and used Spitzer 
to identify stars with 8 /im infrared excesses. In lUzpen et al.l ( 120071 ). we identified 230 stars 
that exhibited 8 /im mid-IR excesses out of 4,015 stars investig ated in the Spitzer Ga lactic 
Legacy Infrared MidPlane Survey E xtraordinaire (GL IMPSE; iBenjamin et al.l |2003| ) and 
MidCourse Space Experiment (MSX: lEgan et al.l 120031 ) catalogs. Out of the 230 stars with 
8 /xm excesses, 183 are of spectral types O or early B. The excess for these stars is plausibly 
explained by free-free emission. Out of the 47 remaining stars of spectral types B8 or later, 
five stars were shown to be fa lse excess sources, and three sources were not yet confirmed 
at 24 yum (jUzpen et al.l 120071 ). In systems with mid-IR excesses we modeled the 8 and 
24 /im excesses as simple blackbodies to characterize the color temperature and fractional 
infrared luminosity of the circumstellar disk. The results were consistent with massive debris 
disks and/or transitional disks with disk temperatures ranging between Ri200-800 K and 
fractional infrared luminosities [Ljn/L^,) ^10~^-10~^. However, infrared photometry alone 
is not sufficient to unambiguously determine the nature of the extraphotospheric emission 
in these sources. Because of their scarcity, identifying new transition disks and warm debris 
disks is an important step toward understanding the evolutionary path of circumstellar disks. 

The nature of the mid-IR excesses emanating from circumstell ar disks around int ermediate- 
mass stars (2-8 Mq) varies greatly and partially depends on age (IZuckermanll200ll ). Young 
stars form circumstellar disks th rough the collap se of the initial molecular cloud as a method 
to conserve angular momentum ( Shu et al.lll987l ). Intermediate-mass pre- main-sequence Her- 
big AeBe stars accre te gaseous material via their circumstel lar disks as they evolve towards 
main-sequence stars ( lHerbig||l960l : IWaters fc Waelkenslll998l ). Circumstell ar dust disks, pre- 
sumably due to planetesima l collisions, surround some main-sequence stars (JBackman fc Paresce 
19931 : iLagrange et al.ll2000l ). Circumstellar gaseous disks that emit excesses owing to free- 
free emission surround other main-sequence stars duri ng late main-sequence or post-main- 



sequence evolution and are known as classical Be stars (IStruvelll93ll : iPorter &: Riviniusll2003 



1999 



McSwain &: Giesll2005l). Circumstellar disks c an surround evolved giants as well (e.g., IJura 



Verhoelst et al.ll2007l : iDeroo et al.ll2007l ). By understanding the nature of the circum- 



stellar disks, we can better understand how these disks evolve. 

The two most common stellar contaminants in surveys for stars with debris disks are 
Herbig AeBe stars and classical Be stars. Herb i g AeBe stars are young stars accreting 



circu mstellar material via their disks (IStromlll972 



Hillenbrand et al. 1992; Hartmann et al. 



19931 ). Emission lines and near-infrared excesses evidence this accretion (JHillenbrand et al. 
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19921 ). IWaters fc Waelkena (119981 ) suggest the following as criteria for Herbig AeBe stars: 
(1) an emission-line A or B star, (2) an infrared excess owing to circumstellar dust, either 
hot or cool, and (3) luminosity class III-V. As the primordial disk dissipates the emission 
lines weaken an d the infrared excess decreases, appearing only at longer wavelengths (e.g., 

y rotating B stars with Ha emission 
1974l:lBalonall2000h . These stars are 



The et al.lll994l ). Classical Be stars are ty pically rapid . 



and an IR excess from free-free emission (IGehrz et al. 



variable in both emission-line strength and luminosity (iDachs et al.lll988l ). The circumstellar 



disks for classical Be stars are thought to be mechanisms for mass loss, unlike Herbig AeBe 
circumstellar disks which facilitate accretion (JBalonall200Cll ). Classical Be stars are sir ailar to 
Herb i g AeBe stars in that they both exhibit emission lines due to circumstellar gas (IStruve 
193ll : lHerbig||l960l ). The greatest distinction between Herbig AeBe stars and classical Be 
stars is the nature of the infrared excess — dust in the former and thermal bremsstrahlung in 
the latter. 

In this paper we present results of optical spectroscopic analysis for 30 mid-IR excess 
sources identified i n the Spitzer GLI MPSE and MSX catalogs to determine the nature 
of infrared excess (jUzpen et al.l 120071 ). By investigating the nature of the infrared excess 
through optical spectral features, we can distinguish between various types of circumstellar 
disks that persist around intermediate-mass stars. Our goal is to remove classical Be and 
Herbig AeBe stars from our sample to identify transition and warm debris disk candidates. 
In §2 we describe how we obtained new optical echelle spectroscopy. We discuss the results 
of our rotational velocity analysis in §3. In §4 we present our methods for determining 
effective temperatures. In §5 we discuss the spectra including the Ha, Fe II A 5317 A, 
Paschen, O I, and Ca II IR triplet line profiles and any implications on the nature of the 
disk from the lines. We describe how we estimate stellar ages through main-sequence fitting 
utilizing our determined fundamental parameters in §6. In §7 we discuss how all of the 
derived parameters provide insight into the nature of the infrared excess, how we use these 
parameters to estimate the contribution of free-free emission associated with our sources, and 
how we draw conclusions regarding the evolutionary state of the circumstellar disk. In §8 
we describe the individual sources and classify the infrared excesses. Finally we use spectral 
classification schemes for stars with infrared excesses, such as classical Be stars, with our 
stars to determine if the classifications are consistent in §9. We also apply our classification 
schemes to well-known stellar samples to determine their classification in §9. 
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2. Target Selection and Observations 



We observed 30 of 42 stellar sources that exhibited an 8 fim mid-lR excess in lUzpen et al. 



(120071 ). 25 of which have an excess at 8 and 20+ /xm. The purpose of these observations 
was to identify circumstellar disk characteristics (e.g., Balmer emission, metallic emission, 
forbidden-line emission), and stellar classification features (e.g. He I, and Mg II absorption) 
based on optical spectral features to determine stellar and circumstellar properties. We chose 
stars with both known and unknown circumstellar evolutionary state. 

We conducted high-resolution, high-signal-to-noise ratio spectroscopic observations at 



the M agellan Telescope with the Magellan Inamori Kyocera Echelle (MIKE; iBernstein et al. 



20031 ) double echelle spectrograph on the nights of 26-28 September 2006 and 7-8 April 2007. 
The spectral coverage was ~3400-9500 A with a resolving power of R=37,000 at 4000 A and 
R=29,000 at 6500 A measured from arc spectra. The typical S/N was ~100:1 in the stellar 
continuum. We conducted standard data reductions in IRAF including fiat fielding with a 
continuum lamp and wavelength calibration using a ThAr arc lamp, but no flux calibration 
was applied. 

We also observed 13 stars to use as rotational velocity standards listed below the solid 
horizontal line in Table [H These stars were used to validate our technique and compare our 
determined values to th ose in the lit e rature. Seven of these stars had projected rotational 
velocities determined in Rover et al.l (l2002ar). Ten of these stars had projected rotational 



velocities determined in ISlettebak et al.l (Il975l ). Four stars had their rotational velocities 
determined in both papers and provide a direct comparison between the two techniques. We 
will discuss these stars further in 53. 



Rotational Velocities 



The most common methods for measuring the projected rotational velociti es are through 



the computation of the fi rst zero of a Fourier Transform of a given line profile (jCarrolllll933 



Grayl Il992l: iRoverl 120051) or thr ough the measurement of a line full width half maximum 
(FWHM) Jsiettebak et al.lll975h . The FWHM of He I A 4471 A, Mg II A 4481 A, and Fe I A 
4476 A are shown to correlate with v sin i in ISlettebak et al.l (119751 ) and is the method we 
adopt to measure projected rotational velocity. However, this method is less reliable at lower 
rotational velocities because of the n on-linearity of the FWH M /projected rotational velocity 
relationship, and the calibration of ISlettebak et al.l (119751 ) h as been shown to b e 10-12% 
lower than measurements made using the Fourier Transform (JRoyer et al.ll2002al ). In order 
to improve upon existing linear relations between He I A 4471 A, Mg II A 4481 A, and Fe I A 
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4476 A we measure the FWHM of the rotationally broadened stellar models of iMunari et al. 
(I2OO5I ) of solar metallicity at a resolving power of R=20,000 from 21,000-7,500 K. 



Figure [T] shows the relationships between projected rotational velocity and the FWHM 
measured over a range of stellar temperatures by the author B.U. of He I A 4471 A (upper 
left), Mg II A 4481 A (upper right), and Fe I A 4476 A (lower left). The He I A 4471 
A line profile includes a Stark broadening component. Since Stark broadening has a non- 
linear dependence on temperature this component degraded the repeatability of consistent 
measurements that caused a greater dispersion in the FWHM measurement. We derived the 
following equations for a linear least squares fit between the measured stellar absorption line 
FWHM and the projected rotational velocity: 



V sin I {km s"^) = 54.25 x He{FWHM) - 46.58 {km s"^), 

V sin I {km s'^) = 47.64 x Mg{FWHM) + 2.01 {km s"^), 
V sin I {km s"^) = 57.04 x Fe{FWHM) - 8.73 {km s"^). 



(2) 
(3) 



The Mg FWHM relation is not linear below ~40 km s~^ due to resolution of the line multi- 
plets or above 300 km s~^ due to blends with the He I A 4471 A line, while the Fe FWHM 
relation is linear at velocities >10 km s~^. 

Table [1] lists the v sin i measurements of our target stars and rotational velocity stan- 
dards and compares them with values listed in the literature, when available. Figure [T]( lower 
right) shows a compar ison b etween our measure ments and literature rotational values from 
Slettebak et al.l (Il975r) and iRover et al.l (l2002al Jbl). Our values are consis tent with those 



derived from 



Royer et al.l (l2002alJbl) and are 10-15% higher than those of ISlettebak et al. 



(I1975I ). This offset compared to ISlettebak et al.l (119751 ) is consistent with results from other 
projected rotational velocity studies. We estimate our uncertainties in projected rotational 
velocity at 10%, given in Table [1], using the FWHM method based on both the dispersion of 
the models from the best-fit line and measurement uncertainties in the FWHM. When two 
rotational velocity indicators are present, the reported rotational velocity is the weighted 
mean of both measurements, and the uncertainty is the uncertainty of the mean. 



4. Surface Temperatures 

We determined effective temperatures by comparing model equivalent width (EW) mea- 
surer nents of He I A 4009 A, He I A 4471 A, Mg II A 4481 A, Fe II A 4233 A, and Ca I A 4227 
A of iMunari et al.l ( 20051 ) R=20,000 solar metallicity models to our spectral measurements. 
For our spectra we used the models in the range 7,000 K < Te// < 21,000 K, with log g'=4.5 
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for Te// < 10,000 K, and log g=4:.0 for Tejf > 10,000 K. We used models with rotational 
broadening between km s~^ and 400 km s~^. The effective temperature step size varies 
from ATe//=1000 K for 12,000 K < Te// < 21,000 K, ATe//=500 K for 10,000 K < T^// 
< 12,000 K and ATe//=250 K for 7,000 K < Te// < 10,000 K. 

The equivalent width ratios of Mg II A 4481 A/He I A 4009 A, Mg II A 4481 A/He I A 
4471 A, Mg II A 4481 A/Fe II A 4233 A, and Mg II A 4481 A/Ca I A 4227 A form monotonic re- 
lations dependent on both rotational velocity and temperature spanning T=21, 000-7,000 K, 
as shown in Figure [2J The solid curves in Figure [2] are polynomials representing the best-fit 
relation between the EW ratio and stellar effective temperature for the models with pro- 
jected rotational velocity of 150 km s~^. The multiple data points at each abscissa/ordinate 
are due to multiple measurements of the same model and show the magnitude of the random 
uncertainties associated with measuring the EW. The dotted curve in Figure [2] (lower left 
panel) represents the best fit polynomial between the EW ratio of Mg II A 4481 A/Fe II A 
4233 A for models with rotational velocity of 250 km s~^ (triangles). We chose Mg and He 
because these lines are prominent in B and early A stars, and He strength decreases rapidly 
in late B stars. Fe lines start to appear in early A stars and increase in strength toward 
later spectral types. Ca lines increase in strength and number starting at late A stars and 
continuing through F stars. At rotational velocities > 150 km s~^ the Mg II A 4481 A/He I 
A 4009 A and Mg II A 4481 A/He I A 4471 A ratios are independent of rotational velocity 
and thus provide the best estimates of Te//. At rotational velocities < 150 km s~^ the Mg II 
A 4481 A/Ca I A 4227 A ratios are independent of rotational velocity and provide the best 
estimates of Tg//. The Mg II A 4481 A/Fe II A 4233 A relationship depends upon rotational 
velocity over all velocity ranges, unlike the other EW ratios. Therefore, this ratio is not gen- 
erally as sensitive to effective temperature, but we use it over a limited temperature range, 
10,000-12,000 K, where other lines ratios are not useful. 

We determined the effective temperature utilizing one or more equivalent width ratios 
in Figure [2] for our target stars and present them in Table [2|, along with spectral types from 
the literature. The uncertainties on effective temperatures were determined by adding in 
quadrature the uncertainty derived from equivalent width fitting with one-half the model step 
size. The former uncertainty component approximately accounts for the random uncertainty 
in measuring the equivalent width, which is dominated by the uncertainty in the definition 
of the stellar continuum. In general, most derived effective temperature measurements are 
consistent with the stellar spectral classification, within uncertainties. 

We found several late B stars to have effective temperatures hotter than their spectral 
classifications would indicate. Such an effect might result if, for example, some stars were 
deficient in heavy metals, resulting in smaller metal-to-helium EW ratios in Figure [2] and 



systematically higher indicated effective temperatures. In order validate the derived tem- 
peratures and test whether this approach is sensitive t o our assumpt i on of solar metallicity, 
we compare in Figures [3] and H] our spectra with the iMunari et al.l (120051 ) stellar models, 
degraded to similar signal-to-noise ratios, in the region near He I A 5876 A. Figures [3] and H] 
show the spectral region near He I A 5876 A (denoted by the dashed vertical line) with the 
best fitting stellar model over-plotted in either red or blue. The narrow strong lines to the 
right of He I A 5876 A are the Na I A 5890/5896 A doublet, generally of interstellar origin 
except in the later type stars which have an atmospheric component. Red curves (24 out of 
30 stars) show the best fit stellar model and indicate that the temperature measured from 
EW ratios fell within one model temperature step size from the best fit model. Blue curves 
(6 out of 30 stars) show the best fit stellar temperature model and indicate that the model 
deviated by more than one model temperature step size from the temperature determined 
from EW ratios. For G014.4239-00.7657 the best fit stellar model was Tej/=18,000 K as it is 
for G299.1677-00.3922 and G007.2369+01.4894. G014.4239-00.7657 and G299. 1677-00.3922 
have the hottest measured Tg// and our relationship for Mg II A 4481 A/He I A 4471 A is not 
very sensitive in this temperature regime. For both stars the Mg II A 4481 A/He I A 4471 A 
derived temperature is ~600 K lower than Mg II A 4481 A/He I A 4009 A ratio pushing the 
weighted mean to a lower temperature. If we ignore the Mg II A 4481 A/He I A 4471 A tem- 
perature, the Mg II A 4481 A/He I A 4009 A temperature, with uncertainties, falls between 
the 17,000 and 18,000 K models, as is found from the best fit models. For the remainder of 
this paper we will use Tejj=16,900±800 K for G014.4239-00.7657 and Te//=16,450±1030 K 
for G299. 1677-00.0392. G007.2369+01.4894 is a known Herbig AeBe star, see §8 for more de- 
tails. The Na I D lines in emission and strong broad He I A 5876 A absorption are indicators 
of circur nstellar gas. This s tar exhibits lin e profile variations consistent with circumstellar 
activity (jCatala et al.lll989l : lPogodinlll994i ). The variable lines and circumstellar gas explain 
the discrepancy between temperatures measured from the He I A 5876 A line fitting and EW 
ratios. G229. 4514+01. 0145 exhibits stronger He I A 5876 A absorption than the Mg II A 
4481 A/He I A 4471 A ratio temperature of Te//=10,220±390 K implies. This is a suspected 
Herbig AeBe star, see §8, and, like G007. 2369+01. 4894, it is better fit with a hotter model, 
Te//=12,000 K. G265. 5536-03. 9951-A and B are binaries and their interactions may effect 
their spectra. For G265. 5536-03. 9951-A we found a better fit with a cooler model of 6,500 K 
rather than the measured T^ff of 7,040+170 K, while G265.5536-03.9951-B is better fit with 
a hotter model of 8,250 K rather than the measured T^// of 7,660+140 K. G307.9784-00.7148 
was also better fit with a cooler stellar model, Te//= 12,000 K rather than the EW ratio mea- 
sured Te//=13, 220+520 K. With the exception of the probable pre-main-sequence stars, the 
best fit models all lie within two stellar models of the measured temperatures based on EW 
ratios. The sources with discrepant model fit temperatures and EW ratio measured tem- 
peratures may show indication of circumstellar gas, see §5. In summary, the six stars that 
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show small discrepancies between temperatures derived from He A 5876 A fits and EW ratios 
are explainable in most cases by circumstellar material and activity rather than metallicity 
enhancements or deficiencies. 

The discrepancy between hterature spectral types and derived temperature for some 
objects in Table [2] may be in part due to differences between observational parameters and 
theoretical stellar models. For example, under the MK system, the definition of a B8 main- 
sequence star requires the He I A 4471 A EW to be equal t o the Mr II A 4481 A EW. Using 
this MK system definition of a B8 main-sequence star and iMunari et al.l (120051 ) models we 
derive a temperature of ^12,600 K through interpolation of the model grid s. Howeve r, the 
observational effective temperature for a B8 main-sequence star is 11,400 K fjCoxll2000l ). We 
do not present revised spectral types for the late B stars in Table[2]due to the small systematic 
temperature difference between optical classification criteria and theoretical models. 



5. Spectral Features 

Metallic and hydrogen emissi on lines in the spectra of both classical Be and Her big AeBe 



stars have been well studied (e.g.. lFinkenzeller &: Mundtlll984l : lAndrillat et al.lll990l ). For ex- 



ample, the lines of He I A 5876 A , Na I D, Ha, Fe II, O I, and Ca II IR triplet have been stud- 
ied a s indicators of circumstella r material around Herbig AeBe stars (iFinkenzeller &: Mundt 



1984; Hamann fc Persson 1992; Boehm fc Catala 1995; Hernandez et al. 



2004h . The spec- 



tra of classical Be stars exhibit, in addition to hydrogen emission, Ca II IR triplet, O I A 



8446 and O I A 7772 A, and Fe II in emission f Andrillat et al. 1990: Jaschek et al. 1993 



Hanuschik et al.lll996l ). By comparing line strengths with stellar models it is possible to 
identify emission and absorption lines arising in circumstellar material and identify the type 
of circumstellar disk. 



5.1. Hydrogen Lines 

5.1.1. Balmer Lines 



The nature of an infrared excess can, in part, be determined through Ha spectroscopy. 
By definition, classical Be stars and Herbig AeBe stars exhibit Balmer emission lines owing 
to their circumstellar gaseous disk. In contrast, debris disks do not exhi bit any Balmer 
emiss ion lines and are devoid of significant amounts of circumstellar gas flLagrange et al. 
2000l ). We measured the Ha equivalent width for our targets and noted the five Ha profile 
types in Table [31 These profile types include single-peaked emission profiles (S), double- 
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peaked emission profiles (D), absorption profiles (A), P-Cygni emission profiles (P), and 
inverse P-Cygni emission profiles (IP). For complex profiles, double-peaked emission, or P- 
Cygni features, we fit multiple Gaussian profiles to the Ha feature and report the sum of the 
EWs in Table El We performed the fitting process manually five times, and we adopt the rms 
dispersi on of the summed EW s as our estimate of the la uncertainty on EW(Hq;). We then 
utilized iMunari et al.l (120051 ) stellar model Ha absorption measurements to correct the Ha 
equivalent width for underlying absorption in order to estimate the true EW of the emission. 
This measure, EW(Hacorr)5 is then used for the remainder of this paper. Normalized Ha 
spectra are shown in Figures I5l&l6l One-third of the sample, ten stars, exhibit Ha absorption 
profiles. Of the remaining twenty stars that exhibit Ha emission, only one star exhibits 
single-peaked emission, one star exhibits inverse P-Cygni, and one star exhibits P-Cygni 
profiles. The majority of the stars that exhibit Balmer emission show a double-peaked 
emission profile (17 stars). This result is not surprising given that the majority of both 
classical Be and Herbig AeBe sta rs exhibit double-peaked Ha profiles (jPorter fc Rivinius 
2OO3I : IPinkenzeller fc Mundtlll984( l. 



5.1.2. Paschen Lines 



Both early- and late-type classical Be stars exhibit Paschen line s in emission with th e 
later type only exhibiting weak emission not reaching the continuum (jAndrillat et al.lll990l ). 
Herbig AeBe stars have been found to exhibit Paschen emission, although the Paschen 
lines are sig nificantly weaker t han the metallic Fe and Ca emission lines in the red spec- 
tral region (ICatala et al.lll986l ). Figures 13 and IH shows spectral regions 8400-8600 A and 
Paschen lines Pal5-19. Lower order Paschen lines either lie off our spectral coverage or 
fall in regions of poor atmospheric transparency. Very few stars exhibit Paschen emission: 
G027.0268+00.7224, G299. 1677-00.3992, G321. 7686+00.4102, and G340.0517+00.6687. Ta- 
ble m lists the spectral features used for confirmation of circumstellar gas and determination 
of the nature of the circumstellar disk. Paschen emission is listed as column 3 in Table HI For 
classical Be stars, emission from metallic lines of Fe II, O I A 8446, O I A 7772 A and Ca II 



is fou nd only when the Paschen lines are in emission (jAndrillat et al.lll990l : iJaschek et al. 



19931). He rbig AeBe stars ma y exhibit O I, Fe II, and Ca II emission with or without Paschen 
emission (ICatala et al.lll986l ). 
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5.2. Metallic Lines 



5.2.1. Call IR Triplet 



For classical Be stars, lAndrillat et al.l ( 19901) found that Ca II is only seen in emis- 
sion if Paschen lines were seen in emission. Ijaschek et al.l ( 19881 ) found that Ca II emis- 
sion was observed in those stars that exhibited a large IRAS excess, >0.6 magnitudes 
at 12 fim. For Herbi g AeBe stars, the Ca II IR triplet is shown to correlate with the 
Ha emission strength (ICatala et al.lll986l ). Only three of our sources exhibited Ca II IR 
triplet emission, G007.2369+01.4894, G027.0268+00.7224, and G229.4514+01.0145, and 
their spectra are shown in Figure [H The best fitting EW ratio determined stellar model 
is shown be neath the stellar spe c trum. G007. 2368+01. 4894 was also shown to exhibit Ca II 



emission m 



Hamann Sz PerssonI (119921 ). and their profile appears similar to our spectrum. 



G229.4514+01.0145 and G027.0268+00.7224 show double peaked Ca II emission. These two 
stars also exhibit double peaked Ha emission. We also find that Ca II emission strength 
is correlated with Ha emission strength, but our sample is not large, containing only three 
sources. Therefore, we do not further investigate this correlation. 



5.2.2. I \ 8446 A Emission and I \ 7772 A Absorption 

O I A 8446 A emission is seen in both Herbig AeBe stars and classical Be stars. The 
O I excitation is due to O I A 1025 A fl uorescence with Lyma. n j3 that subsequently excites 
the related O I A 8446 A ( )Bowenlll947n . kndrillat et al.l (ll99o[ ) found t hat classical Be stars 
only exhibit O I A 8446 A emission if Paschen emission is also present. iJaschek et al.l (119931 ) 
found a correlation between O I A 8446 A emission and O I A 7772 A emission with the 



former about four times stronger than the later, jjaschek et al.l (119931 ) also found that Be 
stars which exhibit both O I A 7772 A in emission also exhibit Fe II A 7712 A in emission. 



Hamann fc PerssonI (119921 ) found that some Herbig AeBe stars exhibit anomalously large 
O I A 7772 A absorption and attribute it to a large circumstell ar envelope. O I A 7772 A 
has also been seen in emission for some Herbig AeBe stars (e.g., iHamann fc PerssonI 1 19921 ). 
Column five of Table H] lists whether O I A 8446 A is found in emission. Column six of 
Table H] indicates whether O I A 7772 A is found in absorption stronger than stellar models 
predict (s) or in emission (e). 



-12- 



5. 2. 3. Fe II Emission 



Classical Be stars exhi bit metallic lines in emission, especially Fe II (jJaschek et al. 



1980l : iHanuschik et al.lll996l ). Presen ce or absence of Fe II is one criterion for classification 
of classical Be stars into subgroups (jJaschek et al.lllQSd ). We show continuum-normalized 
spectra near Fe II A 5317 A in Figures [TUl and [TTl Only four sources exhibit Fe II in emission: 
G027.0268+00.7224 (AOV ) G229.4514+01.014 5 (B9V), G321. 7868+00.4102 (B8VN), and 
G340.0517+00.6682 (B8). Ijaschek et aP Jl98oh Group I sources, stars with Fe II emission, 
are of spectral type B0-B6. These four stars may be an extension to later spectral types 
for Group I classical Be stars. Herbig AeBe stars have also been shown to exhibit Fe I and 
Fe II emission, but the E W of Ha for such sources show strong emission, EW(Ha) < -25 
A jHernandezet al.ll2004h . G340.0517+00.6682 and G321. 7868+00.4102 are the only two of 
the four sources to exhibit strong Ha emission, EW(Hacorr)=-39.38 and -26.81 A and Fe II 
emission. 



5.2.4- Broad Absorption 



Herbig AeBe stars have been shown to exhibit absorption lines such as He I A 5876 
A and O I A 7772 A that a re stronger and broader than expected from the stellar pho- 
tosphere JCatalaet al.lll986h . G007.2369+01.4894 exhibits strong He I A 5876 A as noted 
in §3 and seen in the top left panel of Figure [3l G229. 4514+01. 0145 may also exhibit 
strong He I A 5876 A circumstellar absorption as indicated by the higher temperature stel- 
lar model needed to fit the line than measured by He I and Mg II EW ratios (see §3). 
Several stars exhibit strong broad absorption features at A ~4400, ~4418, and ~4455 A. 
These lines do not appear to correspond with know n diffuse interstellar band s (DIBs) or 
expected photospheric atomic absorption transitions (jJenniskens fc Desertlll994l ). The stars 
that exhibit broad non-stellar absorption lines at these wavelengths are G007. 2369+01. 4894, 
G025.6122+01.3020, G051. 6491-00.1182, and G229.4514+01.0145. These absorption lines 
may imply that G025. 6122+01. 3020 and G051. 6491-00. 1182 are surrounded by circumstellar 
gaseous material similar to the material surrounding Herbig AeBe stars. We will discuss the 
implications of the stellar absorption/emission lines in context of the nature of the infrared 



excess m 
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6. Age Determination 

When stars are not members of known stellar clusters with well-determined ages, stel- 
lar ages are generally estimated by fitting theoretical isochrones to effective te mperatures 



and bolometric luminosities, provided that some distance estimate is available (ISong et al. 



200ll ). Using our derived Teff smd trigonometric pa rallaxes from Hipp arcos /Tycho, we es- 
timated stellar ages by fitting stellar parameters to ISiess et al.l (120001 ) evolutionary model 
tracks. Hipparcos parallaxes were the preferred distance measures and used when available 
(11 stars). In only one instance we used a Tycho parallax. If neither parallax was available, 
we adopted absolute visual magnitudes derived in the literature (4 stars) u sing methods 
such as H/3 luminosity and intri nsic Stromgren colors of luminosity classes (see ISeidensticker 
19891 : IWestinlll985l : lKozoklll985l ). We list parallactic distances in Table [5] except in the cases 
where noted. We then interp olated a bo lometric correction using our derived effective tem- 
peratures in conjunction with ICoxl (J2000l ) effective temperatures and bolometric corrections. 
We applied these bolometric corrections to the Tycho-2 visual magnitudes and correct for 
distance and extinction to derive the absolute bolometric magnitudes. 



Figure [T2] shows our derived absolute bolometric magnitudes and e ffective temperatures 
(i.e., the HR diagram) with isochrones derived from ISiess et al.l (J2000l ) Z=0.02 stellar evo- 
lutionary models over-plotted. Figure [T2] shows the range of stellar temperatures consistent 
with ~2-5 Mq stars. Stars with absolute bolometric magnitudes determined from Hipparcos 
parallaxes are denoted by the diamonds, those with bolometric magnitudes determined from 
Tycho parallaxes are shown by triangles, and those for which literature absolute magnitudes 
were used are shown by asterisks. The zero-age main-sequence is shown by the solid line. 
Isochrones for 1, 10, and 100 Myr are shown by the dotted, dashed, and dash-dot relations. 
The majority of our stars are consistent with 3-5 Mq stellar models. Most of our stars lie 
above the main-sequence indicating either a late pre-main-sequence (few Myr) or late main- 
sequence stellar age (^ 100 Myr) . We estimate a range of possible stellar ages by taking 



the full range of ISiess et al.l (120001 ) model ages consistent with the derived temperatures and 
luminosities, taking into account the measurement uncertainties. Table O lists the possible 
age ranges for the stars with sufficient measurements and compares them to other results in 
the literature. The age ranges we determined for the three stars having published age esti- 
mates were consistent with those measurements. We were able to determine the age ranges 
of two stars identified in literature as pre-main-sequence stars and found that the ages were 
consistent with such a state. Most published measurements of stellar ages have uncertainties 
of order a few % of the stellar lifetime because they are either nearby (small uncertainties in 
luminosity) or derived from main-sequence cluster fitting. Our age measurements typically 
have errors of 10% or greater, so the derived values should be regarded as general age ranges. 
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Implication of Fundamental Parameters on the Nature of the Mid-IR 

Excesses 



Stellar parameters, such as rapid rotation and Ha emission can provide insight into the 
nature of the infrared excess. On average, classical Be stars rotate more rapidly than Herbig 
AeBe stars. Ha emission can be used to determine the amount of free-free emission expected 
from the circumstellar disk. Ha emissi on is also found t o decrease as pre-main-sequence 
stars evolve toward the main-sequence (JManoj et al.ll2006l ). Broadband in frared colors are 



comm only used to discriminate between various circumstellar disks (i.e., iHartmann et al. 
20051 ) . We use the fundamental parameters derived above to determine the nature of the 
circumstellar excess in the following subsections. 



7.1. Large Rotational Velocities 

Intermediate-mass stars rotate at a wide variety of velocities depending on how they 
formed and their evolutionary state. iFinkenzellen (119851) fou n d tha t Herbig AeBe stars rotate 



at values of 100 < v sin i < 225 km s ^. iBoehm fc Catalal (119951 ) predict that higher mass 



Herbig AeBe st ars (> 4 Mp^) shou l d gain angular momentum as they approach the zero-age 



mam-sequence. iBoehm fc Catalal (119951 ) also predict that zero-age main-sequence stars in 
the mass range 2.6 M© < M* < 4.0 Mq should rotate at 205±55 km s~^ and stars with M* 
> 4.0 Mq should rotate at 180±90 km s~^ based on pre-main-sequence projected rotational 
velocities. Classical Be stars are typically rapid rotators having v sin i > 200 km ~^. In 
Figure [12] we compare the rotational velocities for our s ample of 2 ( ) GLI MPSE late-B and 
early A stars (dashed hist ogram), to B8-B9 . 5 star s from lAbt et al.l (120021 ) (solid histogram) 
and classical Be stars from lHanuschik et al.l (119961 ) (dash-dot histogram). Our 20 GLIMPSE 
stars rot ate between 25 km s~^ < v sin i < 310 km s~^ with a mean and dispersion of 209±76 
km s~^. JAbt et al.l (120021 ) found that the rotational velocity distribution for B8-9.5 stars 
is bimodal with peaks at bot h 50 km s~ ^ and 250 km s~^. The mean projected rotational 
velocity and dispersion of the lAbt et al.l (120021 ) late B star sample is 148ib84 km s~^. The 
mean projected rotational velocity and dispersion of iHanuschik et al.l (119961 ) classical Be 
stars is 217±17 km s~^. Thus, our sample has comparatively high rotat ional veloc i ties, m ore 
consistent with the HAeBe stars and Be stars than the field B stars of lAbt et al.l (120021 ). 



We compared our sample rotational velocities to both the classical Be sample and the 
B8-B9.5 sample using a K-S test to determine if our sample could be drawn from either 
of these other distributions. There is no large sample of Herbig AeBe stars with projected 
rotational velocity measurements as they approach the zero-age main- sequence, so we cannot 
make a comparison to these objects. The probability that our sample could be derived from 
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Abt et al.l (120021) sampl e is < 0.1%. The probability that our sample may be derived from 



Hanuschik et al.l (119961 ) sample is 8%. Although th e mean projected rotati onal velocity 



the 

the 

of our GLIMPSE stars is consistent with those of the iHanuschik et al.l (119961 ) classical Be 

stars, it is unlikely that our sample consists of solely Be stars. This result may imply that 

our sample is a mixed group of sources including both normal B stars, classical Be stars, 

and even Herbig AeBe stars. 



7.2. Infrared Colors as a Discriminant for Terrestrial Temperature Dust 



Using the working definition of IWaters fc Waelkend (119981 ) for Herbig AeBe stars, the 
only differing criterion between classical Be and Herbig AeBe stars is the nature of the 
infrared excess. While Herbig AeBe stars have an infrared excess due to circumstellar dust 
and gas, classical Be stars owe their infrared excess to free-free emission in ionized gas. 

Color-color plots are often used to distinguish various types of infrared sources (e.g.. 
Walker et al.lll989l : IZhang et al.l 120041 ). Using a color-color plot in both the near- and mid- 
infrared we can distinguish between sources that have an excess owing to free-free emission 
(and free-bound emission) and sources that have dust emission. In order to model the 
expected cont ribution of fre e- free emission to the infrared continuum we utilize the following 
equation from iTuckerl (119751 ) to estimate the free- free flux. The free-free emissivity at a given 
frequency, v, is 



j^f[y) = 5.44 X 10" 



-39 



QffZ^N^N, 



-hv/kT 



T 



{erg cm ^ s' 



sr 



' Hz-') 



(4) 



where Ne and Np are the electron and proton densities, T is the electron temperature, and 
gff is the Gaunt factor. We assumed an electron temperature of 10,000 K, a pure hydrogen 
gas 



Z=l, and Np=N„ 



=3x10" cm-3. 



The values of electron densities for classical Be stars 

3^ 



vary and we use the mean value from lGehrz et al.l (1l974j ) 



{ue = 3 X 10 cm~ ). Even though 
classical Be stars exhibit such large electron densities within their circumst ellar disks , it ha s 
been shown that such disks are still opt ically thin at optical wa velengths (JMcDavidI l200ll ) . 
We used Gaunt factors from Figure 5 of iKarzas fc Latterl (Il96ll ) by fitting a polynomial to 
parametrize the Gaunt factors as a function of wavelength for T~10,000 K. Equation 5 gives 
an expression for the Gaunt factor for 10,000 K where A is in /^m 



gff = 33.513 - 4.041 x Log{X) + 0.106 x Log{Xf + 0.001 x Log{\f. 



(5) 



In the near-infrared, the contribution of free-bound emission is not negligible. lAshok et al 
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( ll984j ) found the contributions from free-free and free-bound to the total nebular emissi on at 
2.2 yum are 63% and 37%, respectively, for the physical conditions of iGehrz et al.l (Il974l ). We 
therefore include free-bound emission in our calculation of nebular emission. The free-bound 
emissivity as a function of frequency is 



Jx 



fb 



1.4 X 10-^^n-^N^NpT-^/^Z^e^^"-^''^^''^{erg cm~^ s'^ sr'^ Hz' 



(6) 



where n is the excitation level of hydrogen (JTuckerlll975l ). We adopted the same values for 



Ne, Np, and T as for the free-free emissivity above. 

We used the above equations for free-free and free-bound emission to model the infrared 
colors of a two-component spectral energy distribution consisting of a stellar atmosphere 
and a thermal bremsstrahlung excess. We simulated the free-free contribution to the overall 
SED by adding increasing amoun ts of fre e -free (and free-bound) emission (normalized at 8 
/im) to a B8V stellar model from iKuruca (119931 ). We varied the free-free contribution from 
0.01% of the stellar blackbody at 8 /im to 100 times the stellar blackbody at 8 /xm. Varying 
the electron temperature of the circumstellar component between 5,000 K and 15,000 K has 
a negligible effect on the broadband IR colors. We also modeled the colors expected for a 
two-component spectral energy distribution consisting of a stellar blackbody and a cooler 
blackbody representative of circumstellar dust at temperatures 100-1100 K. We then used 
these colors to determine if we can distinguish between free-free emission and dust emission 
in color space. 

Figure [T3] shows a plot of 8-24 //m vs. K-8 /xm colors for the 19 GLIMPSE stars with 
[24] excesses. Open circles represent GLIMPSE sources having Ha in emission and open 
stars represen t GLIMPSE s ource s having Ha in absorption. Classical Be stars (asterisks) 
are those from lZhang et al.l (120041 ). We ut ilized MSX A band and I RAS 25 /xm as 8 /im and 
24 tiui measuremen t s for t he so urces from lZhang et al.l (120041 . |2005| ) respectively. Thirty-five 
of the lZhang et al. ( 2004. 120051 ) stars have MSX A band and IR AS 25 Atm measur e ments . 
Both Be catalogs. IZhang et al.l (|2004|. |2005|). dr a w on Be stars from lJaschek fc Egret! (Il982l ). 
however at least 7 of the 35 iJaschek fc Egret! (119821 ) sources have some indication in the 
literature that they may be Herbig Be stars rather than classical Be stars. The arrow 
shows the reddening vector for Ay=5.0. The solid curve denotes a B8V main-sequence star 
with Av=1.0 and increasing contribution from a free-free component at 8 /im. The solid 
symbols denote where the circumstellar contribution at 8 /im constitutes 5% [diamond), 
33% (hexagon) and 50% (square) of the total 8 /xm flux. The dashed curves denote a B8V 
main-sequence star with Av=1.0 and an 800 K or 300 K blackbody circumstellar component. 
We added the blackbody contribution in a similar manner as the free-free contribution, and 
the solid symbols denote 5%, 33%, and 50% fractional circumstellar contributions. 
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Figure [H] shows that this color-space can distinguish between free-free and cool dust 
origins of the IR excesses. Points to the left of the solid curve may be explained by either hot 
dust or free-free emission. Thus, hot dust is indistinguishable from free-free emission in this 
color space. Points to the right of the solid curve can only be explained by a cool (terrestrial 
or cooler) dust component. The majority (13 of 15) of sources with Ha emission {open circles) 
lie to the left of the solid line, consistent with them being classical Be stars or stars with hot 
circumstellar dust. The majority of the sources with Ha absorption {open stars) lie to the 
right of the solid line, consistent with a cool circumstellar dust component. The majority 
of the classical Be stars {asterisks) lie to the left of the solid line, as expe cted for free-free 



excesses and consistent with their classification. There are, however, 14 iJaschek &: Egret 



([19821) sources that fall to the right of the solid curve. They are HD 6343, HD 37318, HD 
251726, HD 52721, HD 53367, HD 90177, HD 141926, BD-11 4667, HD 174571, HD 204722, 
HD 240010, HD 37967, HD 259431, and HD 50083. Th is sub-sample of 14 star s contains five 



of the seven possible Herbig Be stars within the 35 star lJaschek fc Egretl (119821 ) sample. The 
two remaining possible Herbig Be stars are HD 50138 and HD 101412. These two sources 
show the greatest K-8 color and fall off the top of the plott ed region and to the lef t of the solid 



curve. HD 50138 is also listed as a possible B[e] star (e.g.. lBjorkman et al.lll998l ). Therefore, 
it is possible that some of the stars classified as classical Be stars may actually be Herbig Be 
stars. The two Ha emission GLIMPSE sources that fall to the right of the solid curve are 
G036.8722+00.6199 and G311. 6185+00. 2469. Ha absorption sources that fall to the right of 
the solid curve are G014.4239-00.7657, G311.0099+00.4156, and G339. 7415-00. 1904. These 
five sources are candidates for objects with terrestrial temperature circumstellar dust and 
are discussed in detail below. 



7.3. Individual Free-free Contributions 

A quantitative analysis of the Ha and mid-IR excess fluxes can help distinguish between 
circumstellar dust and ionized gas better than IR colors alone. Specifically, we want to 
determine if the free-free excess in the mid-IR predicted by the Ha flux can account for the 
observed IR excess or if a dust component is required. We use the Ha flux of each star to 
predict the gaseous circumstellar contribution to the 8 /im and 24 /im continuum, under the 
assumption that the Ha, free-free, and free-bound emitting volumes are coincident. Using 
this method one can then verify that the color-color plot in Figure UM discriminates between 
infrared excesses with free- free emission and th ose with an excess due to cool dust. We adopt 



the emissivity at H/3 from iBrocklehurstI (Il97ll ) for 10,000 K, 
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47rJHf3 = 1.24 X IQ-^^NeNp {ergs cm-^ s'^), 



(7) 



where Ne and Np are the number densities of electrons and protons in cm ^ , respectively. 
We used the ratio of Ha to H/3 emission coefficients from iBrocklehurstI (Il97ll ) for Case B, 



jHa 
JHf3 



2.87. 



(8) 



Finally the flux due to free-free emission can be derived utilizing the flux at Ha and the 
emissivity at Ha, 



Fffiu) 



Ha- 



3ff{^) 



jHc 



(9) 



where jff{i^) is the thermal bremsstrahlung emissivity at a given frequency, and JHa is the 
emissivity of Ha. We used an electron temperature of 10,000 K but the exact choice of Ne, 
Np, and T is not critical to the ratio oi jfj/jna- The ratio is independent of gas density and 
only slowly varies with T. 

We calculate the Ha flux by multiplying the absorption-corrected equivalent width of the 
Ha fine, EW(Hacorr), by the reddening-corrected flux at R band. D er eddened R fluxes ar e 
listed in Table [3] computed from R-ba nd magnitud e s of iMonetl (Il998l ) fc lMonet et al.l (120031 ). 
We adopt extinction estimates from lUzpen et al.l (120071 ). which may include circumstellar 
absorption and may be larger than true interstellar extinction. The resultant Ha flux thus 
predicts a maximum possible free-free flux, ignoring possible variability and optical depth 
effects. 

As a r neans of verifying our ca lculations, we used the o bserve d Ha fluxes of H II re- 
gions from Iviallefond etaP Jl983h . Iviallefond fc GossI Jl986h . and Iviallefond eraP Jl986h 
to predict the free-free radio fluxes. We compared our predicted radio fluxes to the measured 
values and found th at the predicted 6 or 21 c m fluxes are a factor of ^2 deflcient relative 
to measured values. Ilsrael fc KennicuttI (Il980l ) found typi cal extinctions of 1-2 mag nitudes 
at visual wavelengths for H II regions. The Ha fluxes from IViallefond fc Gosd (119861 ) do not 
take into account extinction, and this may be part of the reason for the deflcit. Non-thermal 
processes in the radio may also contribute to the radio flux. 

Figure [TBI shows a comparison of the data and models for one star, G036. 8722-00. 4112. 
The modeled free-free component plus stellar blackbody (dash-dot-dot-dot curve) is not 
sufficient to explain the measured excess at 24 /im for this source. A stellar blackbody with 
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additional blackbody due to dust (thin solid curve) fits all the measured data and is the 
likely explanation for the excess. 

We applied our free-free model to all stars with Ha em ission and measured [24] fiux. 



Table [6] gives the measured [24] fiux (from lUzpen et al.l 120071 ) for all 20 GLIMPSE stars and 



derived fiux contribution at [24] from free-free emission for the 15 stars with emission and 
[24] excess. The excesses from six sources are attributable to circumstellar gas since the 
removal of the free-free component results in < 3 cr excesses. We designate such objects as 
Type I sources. The remaining nine stars with emission and [24] measurements have excesses 
> 3 cr after the free-free components were removed. We designate these as Type II sources. 
We note that in none of our sources does the removal of the free-free component reduce the 
8 /im or 24 /im fiux below photospheric levels. In all cases, the correction for circumstellar 
free-free contributions either reduces the 8 and 24 /im measurements to the photospheric 
levels or leaves a significant excess at these bands. 



7.4. Variability 



Both classical Be an d Herbig AeBe stars exhibit variability at Ha (JDachs et al.lll988l : 



Waters fc Waelkenslll998l ). Our sources may also be variable at Ha, and that may explain 
some of the Type II sources that exhibit an excess > 3.0 o at [8.0] and/or [24] when the 
free-free component is removed. In order to assess Ha variability among this sample, we 
investigated Ha profile variations for a sub-sample of four northern Type I or II stars and 
a Balmer absorption star, GOll. 2691+00.4208, G036.8722-00.4112 and G047.3677+00.6199, 
G047.4523+00.5132, and G051. 6491-00. 1182. Ideally all Type II sources would be investi- 
gated for variability but five of the Type II stars are in the southern hemisphere and not ob- 
servable from the Wyoming Infrared Observatory (WIRO) (G300.0992-00.0627, G305.4232- 
00.8229, G307.9784-00.7148, G314.3136-00.6977, and G340.0517+00.6687). 

We obtained optical spectroscopy of all stars with the WIRO-Spec instrument on WIRO 
on the nights of 2005 August 14 and 15 using the 1435 I mm~^ volume-phase holographic 
grating in first order. This instrument is an integral field unit spectrograph consisting of 
293 densely packed fibers in a 19 x 20 configuration with each fiber projecting to 1" on 
the sky. This grating provides a dispersion of 0.7 A pixel"^, a spectral resolution of 1.4 A, 
and wavelength coverage of ~1400 A from ~5750-7150 A. Standard reduction techniques 
were implemented, including fiat fielding using a continuum lamp and wavelength calibration 
using a CuAr arc lamp, but no fiux calibration was performed. We achieved a typical signal- 
to-noise ratio of 100:1. 



-20- 



We also obtained optical spectroscopy using the WIRO-Longslit spectrograph on the 
nights of 2006 September 13 and 14. The WIRO-Longslit is a low-resolution spectrograph 
with a dispersion of 1.1 A pixel"^, a resolution of 2.6 A using a 2" slit width, and covers 
~4400 — 6700 A when using the 600 / mni"^ grating in second order. We used a GG455 
blocking filter to remove third order contamination. We performed the same reduction 
procedures for the longslit data as WIRO-Spec data, and the signal-to-noise ratios were 
similar. 

GOll. 2691+00.4208, G036.8722-00.4112, G051.6491-00.1182 displayed a double-peaked 
emission profile all three nights. G047.4523-I-00.5132 exhibited absorption all three nights. 
G047. 3677+00. 6199 displayed a double-peaked emission profile on the nights of 2005 August 
15 and 2006 September 27. However, on the night of 2006 September 13 the star exhibited 
a single-peaked emission profile. The EW(Hacorr) changed ~25% and may be due to Ha 
variability. A variation in Ha may result in a lower-than-expected calculation of the free-free 
emission. We estimated at what level the Ha EW must vary to account for the excess at [24]. 
We find that the EW(Hacorr) must double for the free-free emission to explain the excess at 
[24]. It is interesting to note that even if the EW(Hacorr) were a factor of two larger, the 
excess is still significant at [8.0] even after removal of the free-free component. 

For the six Type II stars in the southern hemisphere, we are unable to verify any Ha 
variability because we have only one observation of these stars. The Ha variability combined 
with the non-contemporaneous nature of the Ha and IR data, may lead us to over- or 
underestimate the circumstellar free-free contribution to the [8.0] and [24] bands. If the 
circumstellar free-free contribution is underestimated, we expect to see a residual IR excess. 
The nine Type II sources may be examples of this effect. However, we find no instances 
of the converse, i.e., sources where the Ha flux over predicts of the circumstellar free- free 
contribution, thereby creating a deficit at [8.0] or [24]. We find that the six southern Type 
II sources need to exhibit Ha equivalent width variation at two or three times the measured 
levels to explain the excesses as free-free emission. Thus, it is unlikely that the Type II 
sources are to be explained as objects owing their excess to free-free emission. 

Optical depth effects at Ha may also lead to underestimates of the free-free emission. 
The Ha optical depth would have to exceed r^l in order to e xplain the excesses of Type II 



sources by free-free emission. As mentioned above, iMcDavidI (J200ll ) shows that a Be stellar 
disk with an electron density of order 10^^ cm~^ is optically thin, with the maximum r < 1, 
at optical wavelengths. 
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7.5. Physical Nature of IR Excess Sources 

Table [7] lists the derived circumstellar parameters for the nine Type 11 stars with the free- 
free emission component removed. Circumstellar d isk temperature and fractional infrared 
luminosity were determined in the same manner as lUzpen et al.l (120071 ) . Briefly, we flt the 
SED as a single temperature blackbody to model the IR excess. Table [7] also lists the derived 
circumstellar parameters for the four non-emission l ine stars, her e after referred to as Type 
111. We list the circumstellar disk parameters from lUzpen et al.l (120071 ) for the flve MSX 
stars with A > 20 /im observations. 

We searched for correlations between derived circumstellar parameters with stellar pa- 
rameters to provide insight into the nature of the excess. Figure [16] shows circumstel- 
lar disk temperature (upper panel) and stellar effective temperature (lower panel) versus 
EW(HacoT-r)- We found a signiflcant inverse-correlation between EW(Hacorr) and circum- 
stellar disk temperature for Type II stars. The correlation coeflicient of -0.81 for the nine 
Type II sources means there is a < 1 % probability of obtaining such a high degree of correla- 
tion by chance. Such a correlation is expected if the Ha emission is arising within a gaseous 
and dusty circumstellar disk undergoing a clearing of the inner disk. That is, the strength 
of Ha decreases and the dust temperature drops as the inner disk becomes anemic. If the 
excess were due to free-free emission, we would expect the strength of Ha equivalent width 
to be correlated with stellar effective temperature. The hotter stars generate more ion izing 
photons when compared to cooler stars and exhibit a larger excess ( IGehrz et al.lll974l ). We 
flnd no such correlation between stellar temperature and Ha equivalent width for Type I, 
Type II or Type HI stars, see Figure | 



Using our age estimates we investigated whether there is a correlation between EW(Hacorr 
age, and L/^/L* in Figure [T71 Omitting the four pre-main-sequence stars from MSX^ we 
used the nine Type II stars, but only six of the nine have age estimates. We plotted the 
younger age for Type II stars, when available. Using the six Type II sources and the one 
Type HI source with an age, we flnd that age is strongly correlated with Ha equivalent width 
with a correlation coeflicient of 0.97 indicating < 1% probability of obtaining such a high 
probability by chance. Using the o l der ag e when available, there is no correlation between 
EW(HacoT-r) and age. iManoj et al.l (120061 ) found that Ha line strength decreases as Herbig 
AeBe stars age. If those stars are part of the evolved tail to Herbig AeBe stars we would 
expect to see such a correlation between EW(Hacorr) and age. 

Using the above fundamental parameters, above correlations, and spectral features we 
can determine the circumstellar disk classification for our entire sample of 30 mid-IR excess 
sources. We indicate our circumstellar disk classification in Table [71 Pre-main-sequence 
Herbig AeBe stars and T-Tauri stars are indicated by H and TT respectively. Debris disks 
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are denoted by D in Table [71 Stars with an excess owing to free-free emission are classical 
Be stars and listed as B. Finally, stars with a transitional disk system are labeled with T. 
Five of the MSX stars do not have a confirmed excess at A > 20 /xm and thus, we do not 
classify their possible circumstellar disk. We find that utilizing Figure [H] as a first order 
discriminate is useful but may not be the only parameter necessary to determine the nature 
of the excess. Each object is discussed in the following section. 



8. Individual Objects 

In the subsections below, we discuss individual sources. When available, we list the HD 
catalog number, the GLIMPSE or MSX ID number, and the literature spectral type. In 
some cases, the spectral type found in the literature is discrepant from our results. Figure [18] 
is a flowchart illustrating the criteria we use to assign an evolutionary state to each source. 
In summary, we require that our candidate debris disks not exhib it any emission lines and 



exhibit i^-[24] colors < 3.3 magnitudes (e.g.. [Hernandez et al.ll2007l ). Transition disk systems 



may not exhibit Paschen emission lines, may exhibit O I A 8446 A emission, may exhibit O I 
A 7772 A in emission or absorption, may not exhibit Fe II emission, and must have an excess 
not entirely explainable by free-free emission. Classical Be stars must have Balmer emission, 
may exhibit Paschen emission, may exhibit Fe II emission, may exhibit O I emission, and 
must have an excess explainable entirely by free-free emission. Herbig AeBe stars may have 
essentially the same optical spectral characteristics as classical Be stars, and are therefore 
difficult to distinguish based on optical spectra alone. The infrared excess in HAeBe stars 
is dominated by dust rather than free-free emission, therefore knowledge of the SED at 
wavelengths longer than ~10 microns is essential for unambiguous classification. However, 
for our sample, HAeBe stars either exhibit O I A 7772 A absorption and Ca II IR emission 
or exhibit O I emission, Paschen emission, and an excess not entirely explainable by free-free 
emission. 



8.1. Candidate Debris Disk Systems 

HD 231033 (G047.4523-F00.5132)-A0V- As one of the few stars in our sample 
that lacks Ha emission, the derived disk temperature (355±8 K) and Lm/L^ (1.2x10^^) 
make this one of the best candidates for being a warm debris disk system. This source falls 
just to the left of the solid line in Figure [TH in the region where excesses due to either 
free-free emission or hot dust are permissible. The disk temperature and fractional infrared 
luminosity are similar to, but slightly higher than, typical debris disks. Since this star does 
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not exhibit any optical emission lines, we classify the circumstellar disk as a massive debris 
disk. This star exhibits stronger O I A 7772 A absorption than the stellar temperature would 
imply. This may indicate the presence of a small amount of circumstellar gas. This star was 
observed multiple times at Ha with both WIRO and Magellan and was not observed to 
exhibit Ha variability. The age derived for this star (log age(yr)=6. 11-6.37 or 7.58-8.19) 
is consistent with a main-sequence or late-stage main-sequence star, which agrees with its 
circumstellar disk classification. The large uncertainty in age is due to the large uncertainties 
associated with the Tycho parallactic measurement (d=909ly55 pc). 

HD 121808 (G311.0099+00.4156)-A3IV- This star lacks optical emission lines of 
any kind, making it another strong debris disk candidate. Like G047. 4523+00. 5132, it 
exhibits stronger O I A 7772 A absorption than the stellar temperature would indicate, sug- 
gesting that it may also be surrounded by a small amount of circumstellar gas. Circumstellar 
disk parameters (TrfjsA;=315]'^3 K and Lir/ L^.=2.7xl0~'^) are consistent with a warm, mas- 
sive debris disk. There is insufficient data in the literature to make a distance estimate, and 
therefore, an age estimate, for this star. 

HD 151017 (G339.7415-00.1904)-A0- This star does not show any spectral fea- 
tures indicative of the presence of circumstellar gas. Derived circumstellar disk parameters 
(Tdjsfc=311±6 K and L//j/L*=1.0xlO~^) are consistent with a warm, massive debris disk. 
There is insufficient data in the literature to make a distance estimate, and therefore, an age 
estimate, for this star. 



8.2. Candidate Transition Disk Systems 

Since transition disks around intermediate-mass stars have only recently been identified, 
we discuss some general properties of our eight candidate transition disk systems as a whole. 
The four transition disk sources with determined disk temperatures above 500 K exhibit 
Ha emission and O I A 8446 A emission. Two of these four sources exhibit O I A 7772 A 
in emission and two exhibit O I A 7772 A absorption consistent with stellar temperatures. 
Three of the four transition disk candidates with disk temperatures below 500 K exhibit Ha 
emission. Two of these four exhibit O I A 7772 A absorption stronger than the best-matched 
stellar model, and two exhibit O I A 7772 A absorption consistent with stellar temperatures. 
Three of the four stars with disk temperatures below 500 K are also in the region of color- 
color space in Figure [H] not explainable by free-free emission. The change in the O I lines 
from emission in systems with hot circumstellar disks to absorption in systems with cooler 
circumstellar disks may be an evolutionary effect as the disk clears. Hotter circumstellar dust 
would be located nearer the star, and the lack of gaseous emission as the inner disk radius 
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becomes larger may be direct evidence of disk clearing. More transitional disk systems need 
to be identified and studied to provide statistical significance to this trend. 

HD 168246 (G014.4239-00.7657)-B8- This star does not exhibit Ha emission and 
has an intermediate value of projected rotational velocity, 165± 10 km s^^. L/^/L^, for 
this star, 1.3x10"^, is higher than typical for debris disks. The derived disk temperature 
of Trfjsfc=191 K is the coolest in our sample and suggests that the circumstellar material is 
located at a relatively large radius of ~75 AU. We find that this star exhibits anomalously 
strong O I A 7772 A absorption compared to the adopted stellar model, suggesting the 
presence of circumstellar gas. Given the lower temperature derived for this disk and the 
large i^-[24] color (5.15 mag), this disk may be optically thick at [24]. This system may 
contain a transitional disk at the end of its clearing phase. We also find that this star lies 
near a strong 70 /im extended source, raising the possibility that the infrared emission at 
[24] is unrelated to the star. We find that the difference between the position of the strong 
[24] point source and the [8.0] point source is 1.1", i.e., small compared to the 6 " resolution 
at [24]. In order to determine the chance of a random coincidence between the two sources 
we used the [24] source density in this field, 54 stars over an area of 0.01 deg^ (0.1x0.1 deg. 
region), to determine that there are 4.2x10"'* sources per arcsecond^. The probability of a 
[24] source randomly falling within 1.1" of any given [8.0] source is 0.2%. Thus, these two 
sources are probably related. Although unlikely, there is a chance that the excess for this 
source may be due to background cirrus. 

(G036. 8722-00. 4112)-B8V- This star exhibits double-peaked Ha emission but no 
other optical emission lines. Disk temperature (T(ijsfc=427 K) and Lm/L^, (7.2xl0~^) are 
similar to, but slightly higher than, debris disks. Given that this star falls in the region of 
8-24 yum vs. K-9> /im color space not explainable by free-free emission, we classify this star 
as having a transitional disk. The possible ages for this star are widely variant, similar to 
either a late-stage pre-main-sequence star or late-stage main-sequence star. Given the nature 
of the excess and the large difference between the possible ages, we adopt the younger age 
range of log age(yr)=6. 11-6.32. 

HP 180398 (G047.3677-t- 00.6199)-B8Ve- This star has been studied as a classical 
Be star (Jaschek &: Jascheklll992l) . Our Magellan Ha profile shows the same shape recorded 



by Jaschek fc Jaschekl ( 19921 ). but we also observe profile variability consistent with the report 



of iHamdy et al.l ( 199ll ). This star exhibits O I A 8446 A emission but not Paschen emission. 
We classify this source as a transition disk. The derived disk temperature of 784^24 K is 
one of the highest in our sample, as is its Ha equivalent width (EW(Hacorr)=-22.47). The 
determined Lju/L^, of 2.8x10"^ is consistent with a massive debris disk or transition disk. 
The age of this star is somewhere between a late-stage pre-main-sequence star and an early- 
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stage main- sequence star (log age(yr)=6.87-7.72), consistent with a transitional disk. This 
star lies in a region of color-color space explainable either by free-free emission or a hot dust 
excess. 

(G300. 0992-00. 0627)-Al- This star was tentatively associated with an X-ray source 
detected by the Uhuru Satellite (lGurskylll972l ). However, the Rosat Faint Source Catalog 
position for the X-ray source is 580±31 " from the GLIMPSE IR position, effectively ruling 
out an association with the IR excess source (IVoges et al.l 120001 ). This star exhibits strong 
double-peaked Ha emission, O I A 8446 and 7772 A emission, without Fe II or Paschen 
emission. We characterize the infrared excess as due to a transition disk. Circumstellar 
disk parameters of Tdisk=526t.ll K and L/^j/L* =2. 1x10"'^ indicate a disk between those of 
Herbig AeBe stars and those of debris disks. This may be an early form of transitional disk. 



HD 114757 (G305.4232-00.8209)-B6/8V(E)- ISeidenstickeii ( 119891 ) derived My= 
-1.89±0.44, Ay=l.ll±0.13 mag, and distance modulus of 10.03±0.46 mag. Ha exhibits a 
double-peaked emission-line profile. This star exhibits a warm circumstellar disk temperature 
(TdJsA:=675 K) and moderate Lj^/L^ (3.8x10^^). Stellar age estimates indicate either an 
early pre-main-sequence star or a middle/late-stage main-sequence star. This source falls 
in a region of color-color space that may be explained by free-free emission or a hot dust 
component. Free-free emission is not sufficient to explain the large excess at [24], suggesting 
the presence of circumstellar dust. This star exhibits both O I A 8446 and A 7772 A in 
emission without Fe II or Paschen lines in emission. We classify this source as a transition 
disk and adopt the younger stellar age of log age (yr)=5.88. 



HD 118094 (G307.9784-00.7148)-B8V(N)- iLevenhaeen fc Leistei] (120041 ) derived 
the fundamental parameters for this star, which they list as a Be star. They found v 
sin i=239±30 km s-\ Te//=13,500±550 K, log ^=3.90±0.15, log age(yr)=8.03±0.04, log 
L0=2.59±O.O8, and M0=4.2±O.2. Our derived values of projected rotational velocity, 
240±20 km s~^, and temperature, 13,220±520 K, are in agreement with their derived val- 
ues. We determine this star to be either a young pre-main-sequence object or a late-stage 
main-se quence star , log ag e(yr)=6.00-6.23 or 7.66-8.03. Even using their derived parameters 
and the ISiess et al.l (120001 ) models we can not rule out the possibility of a ~1 Myr star. The 



Ha emission profile is highly asymmetric and double-peaked. This star exhibits O I A 7772 
A in emission without Fe II or Paschen lines in emission. Based on our stellar parameters, 
spectral features, Tdisk=557 K, Liii/L^=9.0xl0~*, and the possibility of a few Myr age, 
which we adopt, we classify this excess owing to a transitional disk. 



HD 122620 (G311.6185+00.2469)-B8/9IV/V-[Westin| (119851 ) derived My =-0.613, 
distance=965 pc, Te//=11,190 K, and log age(yr)=8.275. Our derived effective temperature 
of 10,540±255 K is lower. Our age estimates are log age(yr)=6.08-6.18 and 8.16-8.18. The 
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older age is consistent with the hterature but inconsistent with the nature of the excess. The 
free-free component of the [24] excess is one of the smallest of the sample. The circumstellar 
disk parameters (Trfjsfc=306 K and L//j/L*=1.2xl0~^) are consistent with a warm massive 
debris disk. The weak double-peaked Ha emission profile is inconsistent with a debris disk 
nature for this source. We classify this source as a transitional circumstellar disk and adopt 
the younger stellar age range. 

HD 126578 (G314.3136-00.6977)-AlIV- This star exhibits the weakest Ha emission 
profile. The circumstellar parameters (Trfjsfc=328 K and L//j/L^,=9.6xl0~'^) indicate a debris 
disk, but the weak Ha emission precludes such classification. O I A 7772 A is seen in 
absorption with an EW greater than the best-fit stellar model would indicate. This star 
falls in a region of color-color space where excesses due to free-free emission or hot dust are 
permissible. This star may be at the end of its transitional disk phase. 



8.3. Probable Be stars (Type I Sources) 

HD 165854 (G011.2691+00.4208)-B9V- We find that the excess for this star can 
be explained by free-free emission which supports the standing classification of th is star as 



a clas sical Be star. The only emission we find is a weak double-peaked Ha p rofile. iHalbedel 



( 19961 ) found the rotational velocity to be 250±10 km s ^, and lYudinl (120011 ) found it to be 
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242±10 km s ^. These values are slightly greater than our derived value of 220±15 km s 

HD 172030 (G027.0268+00.7224)-A0V- This star exhibits a double-peaked Ha 
emission profile and Fe II emission. The spectrum for this star also shows the Ca II IR 
triplet, O I A 8446 A, O I A 7772 A , and the Paschen lines in emission. The excess for this 
source can be fully explained by free-free emission. It is interesting to note that this star has 
a low projected rotational velocity, 25±5 km s~^, which is not typical for classical Be stars. 

HD 183035 (G051.6491-00.1182)-A0V- The core of the marginally double-peaked 
Ha profile is nearly filled in, and no other emission lines are present. This star does exhibit 
broad absorption features, i.e., A ~4455 A inconsistent with the measured stellar tempera- 
ture. This absorption may be an indication of circumstellar gas. However, this star falls in 
a region of color-color space explainable by free-free emission and has an excess consistent 
with such an origin. We conclude that this is a classical Be star. 

(G299. 1677-00. 3922)-B8- This star exhibits double-peaked Ha emission, Paschen 
emission, and O I A 8 446 A emission. This type of spectrum is typical for classical Be stars 



(lAndrillat et al.lll990l ). 
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HD 107609 (G299.7090-00.9704)-B8/9IV- This star exhibits a large fraction of 
the excess at [24] as free-free emission. However, this source does not exhibit any hnes other 
than Ha in emission. Nevertheless, it is likely that this is a classical Be star. 



HD 121195 (G310.5420+00.4120)-B8IV(N)- IWestinI (Il985f ) derived My=-2.953, 
distance=l,945 pc, and Te//=12,940 K for this star. Their temperature is in agreement 
with our der ived effective t emperature of 12,430±510 K. We determined log age(yr)=~5.30 
utilizing the IWestiru (119851 ) absolute visual magnitude. This extreme youth is probably not 
reahstic given that, at such an age, this star would be embedded and undergo substantial 
extinction, which is not observed. The Ha emission profile is weak and double-peaked, while 
O I A 7772 A is found in absorption. Since the excess can be explained fully by free-free 
emission we classify this source as a Be star, but the optical spectrum is similar to some 
transition disks. 



8.4. Other Circumstellar Disks and Infrared Excesses 



HD 163296 (G007.2369+01.4894)-A1V- This MSX star has been studied numer- 
ous times for both stellar parameters and circumstellar disk evolution. This Herbig Ae (HAe) 
star is o nly 122tjq p c dista nt. Rotational velocities for this source are con sistent in t he lit- 



erature, iFinkenzellerl (Il985l ) derived a rotati onal velocity of 12 0_.^n km s~ , iHalbedell (Il996l ) 
derived a rotational velocity of 120 km s~^, iMora et al.l (1200 ll ) derived a velocity of 130±6 
km s~^, all consistent with our calculated value of 145±15 km s~^. Stellar fundamental 
parameters were calculated to be Te/f=9,200±270 K, MRn;= 0. 48=1=0. 10 from ICidale et al. 
J200lh . and Teff=9,475K, M,=2.4 M©, and 5 Myr old from [Mannings fc Sareentl Jl997h . 



The age of this star varies in the literature but is 4-6 Myr with Ha equ ivalent width vary- 
ing between -12 - -17 A ( Kessler-Silacci et al.l 120061 ). IPogodinI (Il994j ) found this star to 
exhibit rapid Ha variability. Our derived Teff of 10,350=b230 K is hotter than literature 
values. We find that the Ha emission line is stronger than the values reported in the liter- 
ature (EW(Ha)=-26.15±0.15 A). The Mg II A 4481 A line is known to vary, and this may 
explain the differen ce between our measured effective temperature and those in literature 
( iCatala et al.lll989l ). Spectral types f or this star vary between A0-A7, confirming a range 
of stellar temperatures for this star (IFinkenzellerl Il985l : iKessler-Silacci et al.ll2006l ). Our 
measured parameters are consistent with the pre-main-sequence nature of this star. 

HD 167905 (G008.3752-03.6697)-F3V- We classify this MSX star as a T-Tauri 
object based on large Ljpt/L^ (2x10^^), weak Ha emission with a P-Cygni profile, low 
effective stellar temperature (Te//=7260 K), and low rotational velocity. 
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HD 161643 (G009.4973+05.2441)-A7V- We do not have enough data to propose 
an evolutionary status for this MSX star. This star is a binary, and this may explain its 
slightly asymmetric Ha absorption profile (JFabricius fc Makarovll2000l ). The secondary star 
is ^3 magnitudes fainter at V in the Hipparcos catalog. Longer wavelength confirmation of 
the excess is needed. Confirmation of an excess having small Ljr/L^, would indicate a debris 
disk. This star does not exhibit any emission lines. 



Abt&Morre 



ll (Il995r) derived a rotational 



Royer et al.l (J2002bl ) found the rotational 



HD 171149 (G025.6122+01.3020)-B9V- 

velocity for this MSX star of 280 km s^^, while 
velocity to be >301 km s~^. Our derived value of 290±20 km s~^ is consistent with both 
prior measurements. We find Ha in absorption. We also find this source to exhibit broad 
absorption at ~4455 A that may be due to circumstellar material. Additional longer wave- 
length observations are needed to classify this star. The Ha profile needs to be studied for 
variability and emission. Longer wavelength infrared/radio measurements are needed to con- 
firm the excess. With additional optical and infrared/radio measurements it will be possible 
to distinguish between an excess due to a classical Be star, debris disk, or false excess. 



HD 182293 (G054.5163+02.4470VK1IV- Although listed as a KllVin lWright et al 



(I2OO3I ). our spectrum (not shown) indicates a spectral type of II-III. We omit further analysis 
of this source due to its evolved nature. 

HD 58647 (G229.4514-H01.0143 VB9V- There is a d ispa rity between measur ed pro- 
jected rotational velocities for t his star. iGradv et al.l ( 19961 ) and iManoj et al.l (120021 ) found 
it to be 280 km s~^, however, iMora et al.l (120011 ) found the projected rotational velocity 
to be 118±4 km s~^. Historically this star has been classified as a Herbig AeBe star, but 



m 



Manoj et al.l (120021 ) they c onclude that th i s star is a classical Be star based on the low 



fractional infrared luminosity. iMonnier et al.l (120051 ) studied the near-IR and disk properties 
of this source. They measured a stellar luminosity of 295±50 Lq, and found a very low 
fractional infrared luminosity when compared to other Herbig AeBe stars. IVink et al.l (120021 ) 
found that Ha line polarization for this star does not agree with that of classical Be stars but 
exhibits similar properties to Herbig AeBe stars in the (Q,U) plane. They also measured the 
Ha equivalent width to be -8.6 A which is consistent with our measurement of -7.82±0.45 A 
assuming they have similar uncertainties. We find this source to exhibit the Ca II IR triplet, 
O I A 8446 A and Fe II in emission. This source also shows strong O I A 7772 A , and ~4455 
A absorption. Based upon the spectral properties of this source we conclude that HD 58647 
is a Herbig AeBe star. Given the low fractional infrared luminosity, this star may be a more 
evolved form of Herbig AeBe star. 

HD 72106 (G25 7.6236-H00.2459 )-A0IV- This is a visual double star with a separa- 
tion of 0.8 arcseconds (jVieira et al.ll2003l ). We have a spectrum of only the primary star, and 
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no emission lines are found. IVieira et al.l (120031 ) conclude that the secondary of this system 
exhibits the infrared excess, the Ha emission, and is an evolved HAeBe star. The mid-IR 
spectrum of this source exhibits simil ar features of solar s ystem comets also indicating the 
evolved nature of a Herbig AeBe star (ISchiitz et al.ll2005l ). IWade et al.l (120051 ) measured the 
magnetic field of both the primary and secondary of this system and found the primary to 
exhibit a strong magnetic field. They argue that the primary may be a progenitor to the 
Ap/Bp stars. 

HP 73461 (G265.5536-03.99 51)-A5VfcA6/8V- This is a known wide visual bi- 
nary (JBrosche fc Sinachopouloslll989l ). We have spectra of both the primary and secondary. 
Neither source exhibits any emission lines. The primary exhibits stronger O I A 7772 A 
absorption that may indicate circumstellar gas. Longer wavelength confirmation of the ex- 
cess is necessary to determine the nature of the circumstellar disk. Assuming the binary 
system is coeval, the stellar parameters constrain the primary to 2.2±0.2 Mq and secondary 
to 2.0±0.2 Mq with an age of 6.81-6.83 Log yr. 



HD 74534 (G269.5873-05.8882)-G0V- 

emission. We do not detect this in our spectrum. 



Cutispoto et al. 



Cutispoto et al. 



mm detect Ca II A K 
(120031 ) derived a mass of 2 
Mq, an age of < 2 Gyr, and a projected rotational velocity of 17 km s^^. They find this star 
to be an evolved giant. We also find this star to be an evolved giant based on our spectrum. 
Our derived projected rotational velocity of 10±5 km s~^ is the lowest determined for our 
sample. We do not derive additional parameters for this source due to its evolved nature. 



HD 135354 (G321.7868+00.4102)-B8VN- Kozokl (I19851 ) derived distances, abso- 
lute visual magnitudes, and reddenings for Be stars, including this star. The absolute visual 
magnitude is -1.18 at a distance of 1.01 kpc. Utilizing these parameters, we determined a 
stellar age range of log age(yr)=5.83-6.12 and 8. 01-8.13. This s t ar exh ibits Ha, O I A 8446 
A, O I A 7772 A, Paschen, and Fe II emission. lAndrillat et al. ( 19901) found t hat c lassical 
Be stars exhibiting O I emission also show Paschen emission. iJaschek et al.l (119931 ) found 
that classical Be stars exhibit both O I and Fe II emission. This star exhibits both of these 
characteristics. The fractional infrared luminosity and disk temperature (L/r/L*=2.1x10~'^ 
and TrfjsA:=556 K) are consistent with transition disks but the measured Ha EW is the second 
strongest of the sample. Given that the Ha flux is insufficient to explain the excess at [24] 
and the EW of Ha is consistent with Herbig AeBe stars we classify this star as a Herbig 
AeBe star. 

HD 151228 (G339.4392-00.7791)-A0/l IV/V- The excess at [8.0] is not confirmed 
for this star at [24] , making it a false mid-IR excess. 



HD 150625 (G340.0517+00.6687)-B8- iKozokl (119851 ) derived an absolute visual 
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magnitude of -1.87 and a distance of 0.83 kpc. Using these parameters we constrain stellar 
ages to log age(yr)=5.62-5.85 and 7.81-7.90. This star exhibits a strong Ha emission profile, 
Fe II, Paschen, and O I A 8446 and 7772 A emission lines. The measured Ha flux is insufficient 
to explain the excess at [24]. This is the other Type II source that we classify as a Herbig 
AeBe star based on the emission line spectrum. The fractional infrared luminosity and 
disk temperature (L/ij/L^,=1.5xl0~^ and Tdisk=823 K) for this star is the largest in the 
GLIMPSE sample and falls between Herbig AeBe stars and debris disks. This may be an 
evolved Herbig AeBe star. 

HD 152404 (G347.3777+04.2010)-F5V- This sta r is a spectroscopic binary in the 
Upper Centarus-Lupus with an orbit al period of 13.6 d ays (JManset et al.ll2005l ). This star is 
~ 17 Myr old and classified as F5V (jChen et al.ll2005l ). This star has a st rong excess ratio 
given its age and may be a long lived primordial disk (jSiegler et al.l 120071 ). Our projected 
rotational velocity fo r the two stars is 15±5 km s"^ consistent with 18.5±1.0 km s~^ from 
(jAlencar et al.ll2003l ). Due to the nature of this binary source we do not derive effective 
temperature since our method of deriving temperature does not extend to spectroscopic 
binaries of similar mass. The Ha profile indicates an unusual double peaked inverse P-Cygni 
profile. This may indicate circumstellar accretion. 



9. Sample Comparisons 

Given the similarity of the Ha profiles for both stars classified as candidate transition 
disks systems and classical Be stars, we make a comparison of the spectral f eature s for 
our sample w i th a g roup of well-studied classical Be stars f rom lAndrillat et al.l (Il990l) and 
Jaschek et al.l (119931 ). We use the following conclusions from lAndrillat et al.l (Il990l ) for clas- 
sical Be stars: 1) Ca II emission is infrequent for spectral types later than B6 and nearly all 
stars that exhibit Ca II in emission exhibit Paschen emission, 2) Paschen lines are in emission 
when O I is in emission, 3) When IRAS 12 yum is in excess greater than 0.55 magnitudes 
of the predicted measurement then Ca II is in emission. In order to make the comparison 
with the last criterion we used the [8.0] measurement in excess greater than 0.55 magni- 
tudes. Table [13] lists the 10 stars we are able to classify as either consistent or inconsistent 
with a classical Be star classification based on the three criteria we mentioned above. In 



all 10 instances the classifications we derive using the criteria of lAndrillat et al.l (Il990l ) are 
consistent with our classifications. 

In order to validate our classi fication scheme we als o applied the classifications in Fig- 
ure [18] to the classical Be stars in lAndrillat et al.l (Il990l ). In order to make the comparison 



over all spectral features classified, we supplemented the spectra from lAndrillat et al 



(119901) 
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with spectra fro m I Jaschek et al. 



19 stars in both Andrillat et al 



( 1993 ) to g ain the O I A 7772 A feature. There were only 
(Il990l ) and IJaschek et al.l (Il993l ). Using our classification 
criteria we were able to classify only one of the 19 stars as a classical Be star based on spectra 
alone. The remaining 18 stars required some knowledge of the SED. We fit Kurucz stellar 
models to 2MASS photometry for these 18 stars in a manner similar to our GLIMPSE stars 
in order to estimate the amount of extraphotospheric flux emanating from these sources. 
Since these 18 stars are all bright (~V< 5) and 2MASS photometry for such stars have 
large uncertainties, we chose not to determine the free-free emission from Ha flux because 
of the larger normalization uncertainties associated with the K band. In order to deter- 
mine the nature of the excess emission we determined the slope of the excess component, 
/?, where Fj, oc v^ after an estimate of the stellar flux was removed. The slope of free-free 
emission in the infrared varies from (3 ~ -0.1-2.0 depending on the optical depth of the gas 
(IWright fc Barlowlll974j ). By contrast, sources with dust components would have slopes (3 < 
-0.1. In o rder to determine t he slope of the infrared excess, we used flux measurements at 8.7 
/zm from iGehrz et al.l (119741 ) and the IRAS Point Source Catalog at 25 /im measurements. 
Only 10 of the 18 stars had photometry at both 8.7 and 25 /im. We found that eight out 
of the 10 stars exhibit positive slopes consistent with optically thick free-free emission, r] 
Tau and 48 Per exhibited negative slopes. Only 48 Per exhibits a slope less than -0.1 at the 
3(7 level. Given the length of time between the two infrared measurements (> 10 years) the 
large negative slope may be a consequence of the star's known variability (e.g.. lKoen fc Eyer 
2OO2I ) . Assuming that slopes consistent with free-free emission would cause the stars to be 
classified as Type I sources (Be stars), 10 of 11 known classical Be stars would be classified as 
Be stars using our classification methods. By comparison, all three of our candidate debris 
disk systems exhibit slopes less than -0.1 at a 3cr or greater level. None of the six classical Be 
stars within our sample exhibit slopes less than -0.1 at a 3o" level, meaning that their excesses 
are consistent with optically thin or thick free-free emission. Neither of the two candidate 
Herbig AeBe stars exhibit a slope less than -0.1 at the 3cr level. Four of the eight candidate 
transition disk systems exhibit a slope less than -0.1 at greater than a 3(j, level indicative of 
dust. These four stars are G014.4239-00.7657, G036.8722-00.4112, G311.6185+00.2469, and 
G314. 3136-00. 6977. The explanation for the four remaining candidate transition disk sys- 
tems not exhibiting a /? < -0.1 at a 3cr level may be variability. An additional scenario for the 
shallowe r slopes is a narrow ring of single temperature hot dust, something analogous to C, 
Leporis (jChen fc Jurall200ll ). Another possible explanation for the excess for these four stars 
may be optically thick free-free emission. In order to be conservative, we added a possible Be 
classification to the remaining four candidate transition disk systems, G047. 3677+00. 6199, 
G300.0992-00.0627, G305.4232-00.8229, and G307.9784-00.7148 in Tabled 
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10. Conclusion 

In an effort to identify the nature of the infrared excesses, we obtained high-resolution 
optical spectroscopy for 31 luminosity class V or IV stars of spectral type B8 or later ex- 
hibiting 8 /xm mid-IR excesses in the GLIMPSE and MSX catalogs. We measured stellar 
effective temperatures for 28 stars and projected rotational velocities and Ha equivalent 
widths for 30 stars. The projected rotational velocities of the 20 GLIMPSE stars lie between 
25 and 310 km s^^ with a mean and dispersion of 209±76 km s^^ consistent with the mixed 
nature of the sample. We used derived stellar parameters and distances to estimate stellar 
ages ranging from ~1 Myr to 100 Myr for 16 stars (9 of the 20 GLIMPSE stars). Twenty of 
the 31 stars (16 of the 20 GLIMPSE stars) exhibit Ha emission, consistent with characteris- 
tics of pre-main-sequence stars (i.e., Herbig AeBe stars) or classical Be stars. Ha equivalent 
widths range between -39.38 A (characteristic of classical Be stars or Herbig AeBe stars) and 
9.41 A (main- sequence B and A stars), suggesting a variety of origins in a sample of objects 
spanning a range of stellar evolutionary states. 

For 16 of the 20 GLIMPSE stars showing Ha emission, we use the Ha flux to estimate 
the expected contribution of (free-free and free-bound) emission to the mid-IR excesses. In 
six out of the 16 cases, the expected IR excess can fully account for the observed infrared 
excess (Type I sources). We did not confirm the infrared excess at 24 /im in one case. In 
the remaining nine cases (Type II sources) the expected contribution from free-free emission 
cannot account for all of the observed IR excess, implying that either 1) the circumstellar 
excess is time- variable and the non-contemporaneous nature of Ha and mid-IR data render 
this estimate invalid, or 2) there is a dust component to the IR excess. The lack of evidence 
for Ha variability for the majority of our sources, combined with the fact that in no instance 
does the Ha over predict the IR excess component, makes the first scenario unlikely. We 
conclude that the second scenario, the presence of a circumstellar dust component, is most 
probable. Combining correlations of spectral features for sources with known characteristics 
we find that the nine Type II sources with IR excesses require a dust component and are 
likely to be young main-sequence or near main-sequence objects harboring warm (300-800 
K) circumstellar dust. The nominal fractional infrared luminosities for these nine sources 
after removal of the free-free component is ~10~^-10~^. For the six of these nine sources 
with age estimates, there is a correlation between estimated age and Ha EW consistent with 
a progression of disk clearing with time. Such characteristics make these stars candidates for 
being "transition disk" or "evolved-HAeBe" systems containing rare dissipating primordial 
gas and dust disks. One of the four systems lacking Ha emission has a very large i^-[24] 
color, > 4 magnitudes. Such a large color difference implies that the excess is optically thick 
at [24]. We conclude that this source is also a transition disk. The remaining three systems 
lacking Ha emission have disk temperatures of 300-400 K and fractional infrared luminosities 
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of 10 ^. These characteristics make them candidates for dusty debris disk systems. 

We show that -ft'-S, 8-24 color-color diagrams are a useful tool to distinguish cool (T^^ < 
300 K) dusty mid-IR excess objects from hot (T^ > 500 K) dusty or free-free mid-IR excesses. 
Such diagrams reveal that additional information beyond near/mid-IR colors is necessary to 
distinguish hot dust from free-free excesses. This diagram also reinforces the idea that 
some stars cataloged as classical Be stars may in fact be late-stage pre-main-sequence stars 
containing circumstellar dust. 
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Table 1. Projected Rotational Velocities 



Spec. 



V sin i Lines Used 



(km s^l) 



V sin i (Literature) References 
(km s~l) 



G007.2369+01.4894 


163296 


AlV 


145±15 


M 


G008. 3752-03. 6697 


167905 


F3V 


35±5 


F 


G009. 4973+05. 2441 


161643 


A7V 


140±15 


M 


GOll. 2691+00. 4208 


165854 


B9V 


220±15 


H,M 


G014.4239-D0.7657 


165854 


B8 


165±10 


H,M 


G025. 6122+01. 3020 


171149 


B9V 


290±20 


H,M 


G027.0268+00.7224 


172030 


AOV 


25±5 


M 


G036. 8722-00. 4112 




B8V 


295±20 


H,M 


G047.3677+00.6199 


180398 


B8Vo 


310±25 


H,M 


G047.4523+00.5132 


231033 


AOV 


150±10 


H,M 


G051. 6491-00. 1182 


183035 


AOV 


270±20 


H,M 


G229. 4514+01. 0145 


58647 


B9V 


295±20 


H,M 


G257.6236+00.2459 


72106 


AOIV 


240±25 


M 


G265. 5536-03. 9951-north(A 


73461 


A5V 


40±5 


F 


G265. 5536-03. 9951-south{B) 




A6/8V 


25±5 


F 


G269. 5873-05. 8882 


74534 


GOV 


10±5 


F 


G299. 1677-00. 3922 




B8 


130±10 


H,M 


G299. 7090-00. 9704 


107609 


B8/9IV 


175±15 


H,M 


G300. 0992-00. 0627 




Al 


255±20 


H,M 


G305. 4232-00. 8229 


114757 


B6/8V{E) 


230±20 


H,M 


G307.9784-00.7148 


118094 


B8V(N) 


240±20 


H,M 


G310. 5420+00. 4120 


121195 


B8IV(N) 


270±20 


H,M 


G311. 0099+00. 4156 


121808 


A3IV 


235±25 


M 


G311. 6185+00. 2469 


122620 


B8/9 IV/V 


265±20 


H,M 


G314.3136-00.6977 


126578 


AlIV 


205±20 


M 


G321. 7868+00. 4102 


135354 


B8VN 


230±20 


H,M 


G339. 4392-00. 7791 


151228 


AO/1 IV/V 


40±5 


M 


G339. 7415-00. 1904 


151017 


AO 


240±15 


H,M 


G340. 0517+00. 6687 


150625 


B8 


230±15 


H,M 


G347.3777+04.2010 


152404 


F5V 


15±5 


F 



250±10, 242±10 (1), (2) 

280, 301 (3), (4) 



118, 280 



(5), (6) 



C7) 



22252 


B8V 


255±20 


M 


249 


(8) 


48757 


AOV 


30±5 


M 


27 


(8) 


50853 


AOV 


205±20 


M 


195 


(8) 


181454 


B8V 


80±10 


H,M 


75 


(9) 


100673 


B9Ve 


150±15 


M 


160,125 


(8), (9) 


97277 


A2IV 


55±5 


M 


45 


(9) 


83953 


B6V 


285+20 


H,M 


260 


(9) 


95370 


A3IV 


95+10 


M 


86, 90 


(8), (9) 


83754 


B5V 


190+20 


H,M 


150 


(9) 


98718 


B5Vn 


325+35 


H 


280 


(9) 


90994 


B6V 


105+10 


H,M 


80 


(9) 


78045 


A3V 


50+ 5 


M 


34, 25 


(8), (9) 


83446 


A5V 


135+15 


M 


133, 125 


(8), (9) 



Note. — See lUzpen et al.l J2007D for further information about GLIMPSE or MSX Catalog ID. Line 
profiles used to derive projected rotational velocities are M for Mg II A 4481, H for H I A 4471 A and 
F for F e I A 4476. (11 iHalbedell Jl996|): (21 [Yu din 12001); (3) Abt fc Morrelll i 19951 '): (4) 

2002d l: f5') lMora et al.l j20oil): ( 6') IGradv eTal] 111998 ): (7) .Cutispoto et al.l i2003ll : (8) 

2002al i: Ol lSlettebak et al.l 119751 1 



Rover et al 



Rover et al 
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Table 2. Effective Temperatures 



ID 


Spec. 


(K) 


(1) 


(2) 


(3) 


G007.2369+01.4894 


AlV 


10350+230 


G008. 3752-03. 6697 


F3V 


7260+160 


G009. 4973+05. 2441 


A7V 


8320+130 


GOll. 2691+00. 4208 


B9V 


12350+520 


G014. 4239-00. 7657 


B8 


16460+640 


G025. 6122+01. 3020 


B9V 


12520+510 


G027.0268+00.7224 


AOV 


10510+300 


G036. 8722-00. 41 12 


B8V 


11740+260 


G047.3677+00.6199 


BSVc 


15700+600 


G047.4523+00.5132 


AOV 


12260+510 


G051. 6491-00. 1182 


AOV 


13100+530 


G229. 4514+01. 0145 


B9V 


10220+390 


G257.6236+00.2459 


AOIV 


11530+300 


G265. 5536-03. 9951-A 


A5V 


7040+170 


G265. 5536-03. 9951-B 


A6/8V 


7660+140 


G269. 5873-05. 8882" 


GOIII 




G299. 1677-00. 3922 


B8 


16170+730 


G299. 7090-00. 9704 


B8/9IV 


12520+505 


G300. 0992-00. 0627 


Al 


13510+510 


G305. 4232-00. 8229 


B6/8V(E) 


13520+510 


G307.9784-00.7148 


B8V (N) 


13220+520 


G310. 5420+00. 4120 


B8IV(N) 


12430+510 


G311. 0099+00. 4156 


A3IV 


9800+130 


G311. 6185+00. 2469 


B8/9 IV/V 10540+255 


G314. 3136-00. 6977 


AlIV 


9870+130 


G321. 7868+00. 4102 


B8V N 


12340+505 


G339. 4392-00. 7791 


AO/1 IV/V 8800+150 


G339. 7415-00. 1904 


AD 


12220+510 


G340. 0517+00. 6687 


BS 


14080+520 


G347.3777+04.2010'' 


F5V 





Note. — (a) This star is a giant. Sec §8 
for further details, (b) This star is a spec- 
troscopic binary. Equivalent width measure- 
ments vary too greatly to determine temper- 
ature. See §8 for further details. 



Table 3. Ha Equivalent Width Measurements 



ID 

(1) 


Spec. Type 
(2) 


Date 
(3) 


EW{Ha) 

(A) 

(4) 


EW{Ha„i,J 

(A) 

(5) 


EW(Hacorr) 

(A) 

(6) 


Ro Flux 

(erg cm^ s^ A^ ) 

(7) 


Ha Flux 

(erg cm— 3— ) 

(8) 


Ho, Profile 
(9) 


G007.2369+01.4894 


AlV 


20060928 


-26.15+0.15 






1.42x10"^"- 






S 


G008. 3752-03. 6697 


F3V 


20060928 


-0.17±0.08 






2.12x10"^^ 






P 


G009. 4973+05. 2441 


A7V 


20060928 


5.11±0.02 






2.52x10-1* 






A 


coil. 2691+00. 4208 


B9V 


20060927 


-4.24±0.16 


9.62 


-13.86+0.16 


1.30x10-1^ 


1.80x10' 


-11 


D 


G014. 4239-00. 7657 


B8 


20060927 


4.45±0.03 




0.00 


1.40x10-13 






A 


G025. 6122+03. 3020 


B9V 


20060928 


6.42+0.09 






5.99x10-1^ 






A 


G027.0268+00.7224 


AOV 


20070407 


-1.97+0.43 


11.89 


-13.86+0.43 


1.74x10"!^ 


2.41x10" 


-11 


D 


G036. 8722-00. 4112 


B8V 


20060927 


-2.39+0.15 


12.68 


-15.07+0.15 


1.68x10-1^ 


2.53x10" 


-11 


D 


G047. 3677+00. 6199 


B9Ve 


20060927 


-14.08+0.24 


8.39 


-22.47+0.24 


1.79x10-1^ 


4.02x10" 


-11 


D 


G047. 4523+00. 5132 


AOV 


20060927 


6.49+0.02 




0.00 


5.12x10-13 






A 


G051. 6491-00. 1182 


AOV 


20060928 


5.76+0.15 


8.65 


-2.89+0.15 


1.30x10-1^ 


3.76x10" 


-12 


D 


G229. 4514+01. 0145 


B9V 


20060926 


-7.82+0.45 






8.66x10-1^ 






D 


G257. 6236+00. 2459 


AOIV 


20060926 


9.41+0.06 






1.03x10-1^ 






A 


G265. 5536-03. 9951-A 


A5V 


20060927 


3.77+0.06 






2.36x10-1^ 






A 


G265. 5536-03. 9951-B 


A6/8V 


20060927 


4.93+0.11 






7.01x10-13 






A 


G269. 5873-05. 8882 


GOIII 


20060926 


2.12+0.03 






1.04x10-1^ 






A 


G299. 1677-00. 3922 


B8 


20070407 


-4.37+0.12 


9.62 


-13.99+0.12 


5.51x10-13 


7.71x10" 


-12 


D 


G299. 7090-00. 9704 


B8/9IV 


20070407 


-1.99+0.16 


8.13 


-10.12+0.16 


5.36x10-13 


5.42x10" 


-12 


D 


G300. 0992-00. 0627 


Al 


20070407 


-13.89+0.32 


7.76 


-21.65+0.32 


2.41x10-13 


5.22x10" 


-12 


D 


G305. 4232-00. 8229 


B6/8V(E) 


20070407 


-14.22+0.42 


7.76 


-21.98+0.42 


7.34x10-13 


1.61x10" 


-11 


D 


G307. 9784-00. 7148 


BSV(N) 


20070407 


-4.63+0.69 


8.64 


-13.27+0.69 


1.62x10-1^ 


2.15x10" 


-11 


D 


G310. 5420+00. 4120 


B8IV(N) 


20070407 


-0.09+0.09 


9.88 


-9.97+0.09 


6.88x10-13 


6.86x10" 


-12 


D 


G311. 0099+00. 4156 


A3IV 


20070407 


5.21+0.03 




0.00 


1.67x10-1^ 






A 


G311. 6185+00. 2469 


B8/9IV/V 


20070407 


1.96+0.23 


12.86 


-10.90+0.23 


3.12x10-13 


3.40x10" 


-12 


D 


G314. 3136-00. 6977 


AlIV 


20070407 


4.20+0.29 


16.83 


-12.63+0.29 


9.41x10-13 


1.19x10" 


-11 


D 


G321. 7686+00. 4102 


B8V N 


20070407 


-16.93+0.53 


9.88 


-26.81+0.53 


5.88x10-13 


1.58x10" 


-11 


D 


G339. 4392-00. 7791 


AO/IIV/V 


20060927 


3.28+0.76 


14.48 


-11.20+0.76 


8.58x10-13 


9.61x10" 


-12 


D 


G339. 7415-00. 1904 


AO 


20060927 


5.50+0.03 




0.00 


4.15x10-13 






A 


G340. 0517+00. 6687 


B8 


20060927 


-31.62+0.08 


7.76 


-39.38+0.03 


1.17x10-1^ 


4.61x10" 


-11 


D 


G347. 3777+04. 2010 


B8 


20060927 


-0.02+0.10 






2.08x10-1^ 






IP 



00 



Note. — (5) Hq: equivalent ■width of the nearest stellar model. This equivalent width was then used to determine EW(HcKcorT-) and subsequently calculate the Ha 
flux. (7) The reddening corrected R flux. (8) Hck flux was determined only for GLIMPSE mid-IR excess stars with 24 /im measurements and Ha emission. (9) D: 
double-peaked profllc, A: absorption only profile, S: single emission only proflle, P: P-Cygni profllc, IP: inverse P-Cygni profile. 
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Table 4. Spectral Features 



ID Ha Emission Paschcn Emission Ca II Emission O I 8446 Emission O I 7772 Absorption Fc II Emission 

(1) (2) (3) (4) (5) (6) (7) 



G007.2369+01.4894 


y 


G008. 3752-03. 6697 


y 


G009. 4973+05. 2441 


n 


GOll. 2691+00. 4208 


y 


G014.4239-00.7657 


n 


G025. 6122+01. 3020 


n 


G027.0268+00.7224 


y 


G036. 8722-00. 4112 


y 


G047.3677+00.6199 


y 


G047.4523+00.5132 


n 


G051. 6491-00. 1182 


y 


G229. 4514+01. 0145 


y 


G257.6236+00.2459 


n 


G265. 5536-03. 9951-A 


n 


G265. 5536-03. 9951-B 


n 


G269. 5873-05. 8882" 


n 


G299. 1677-00.3922 


y 


G299. 7090-00. 9704 


y 


G300. 0992-00. 0627 


y 


G305. 4232-00. 8229 


y 


G307.9784-00.7148 


y 


G310. 5420+00. 4120 


y 


G311. 0099+00. 4156 


n 


G311. 6185+00. 2469 


y 


G314.3136-00.6977 


y 


G321. 7868+00. 4102 


y 


G339. 4392-00. 7791 


y 


G339. 7415-00. 1904 


n 


G340. 0517+00. 6687 


y 


G347.3777+04.2010'' 


y 



N/A 



N/A 



Note. — (2) Stars that exhibit Ha in emission are denoted by (y). (3) Stars with Paschen lines in emission are denoted by 
(y). (4) Stars that exhibit Ca II IR triplet emission are denoted by (y). (5) Stars that exhibit O I A 8446 A in emission are noted 
by (y) (6) Stars that are shown to have O I A 7772 A in absorption stronger than stellar models predict are denoted by (s), those 
that exhibit it in emission (e), and those that exhibit absorption are denoted by (n). G269. 5873-05. 8882 and G347. 3777+04.2010 
are listed as N/A because we do not determine their stellar temperature and cannot make a comparison to stellar models. (7) 
Stars that exhibit Fe II emission are denoted by (y). (a) Spectral features indicate that this star is a giant. See §8 for further 
details, (b) This star is a spectroscopic binary. See §8 for further details. 
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Table 5. Estimated Stellar Ages 





ID 

(1) 


Distance 

(pc) 
(2) 


Adopted Mg^i 

[mag] 

(3) 


Age Derived 

(log yr) 

(4) 


Age in Literature 

(log yr) 

(5) 



G007. 2369+01. 4894 
G009. 4973+05. 2441 
G025. 6122+01. 3020 
G036. 8722-00. 4112 
G047. 3677+00. 6199 
G047.4523+00.5132 
G051. 6491-00. 1182 
G229. 4514+01. 0145 
G257. 6236+00. 2459 
G265. 5536-03. 9951-A'' 
G265. 5536-03. 9951-B'' 
G305. 4232-00. 8229= 
G307. 9784-00. 7148 
G311. 6185+00. 2469" 
G321. 7868+00. 4102" 
G340. 0517+00. 6687" 
G347. 3777+04. 2010 



,+ 17 
■^-13 
+32 
■22 

., + 12 



137 
1142 

^"^^-142 

353_gg 

909;?^, 

345^17 

1391^9 
1391^9 

„,„ + 2272 

"^*'-326 

965 

1010 

830 



-l.W 



-0.10'* 



^+1.34 
" -0.34 

-'■^'-0.45 
10.34 
-0.36 
, t.o + 3.40 
-l.oo_j gg 

,;+0.78 

" -1.06 

4-4.01 

-0.37 

10.79 

-1.03 

- + 0.59 

" -0.50 

,+0.84 

-1.26 

r,+0.42 

-0.49 

,; + 0.42 

-0.50 



-2.05 
-2.00 
-0.98 
-2.25 
-0.13 

0.70 

1.85 
-3.05+0.45 
-3.04 



1.38 
2.10 
1.13 ±0.50 
■2.09+0.51 
■ 3.14+0.14 



6.13- 
6.65- 



6.09 
6.26- 



5. 
6.00. 



5.83 
5.62. 



6.63, 
6.95, 

7.09- 
6.32, 

6.87- 
6.37, 

6.13- 
6.43, 
6.55, 

6.81- 

6.81- 
J8, 7, 
6.23, 
6.18, 
6.12, 
5.85, 



8.42 
8.09 
8.10. 
7.72 
7.58 
-8.05 
8.37 
7.33 
-6.83 
-6.83 
87-7. 
7.66 
8.16 
8.01 
7.81. 



■8.27 
■8.74 



8.20 
8.19 



8.43 
8.24 



8.03 8.03+0.04" 
8.18 8.275-f 

8.13 
7.90 



Note. — (a) jMannings 8z Sa rgent 199^); (b) Assuming botli stars to be coeval 
only a 2.2 and 2.0 Mc^ b inary is possible given temperature measurements, (c) Uti- 
lizing ISeidensticken J 19891) extinction and absolute visual magnitude measurements (d) 
jLc venhagcn & Leister 20 0^ (e) U tilizinglWestin ( 1985) absolute visual magnitude, (f) 
IWestin, 1 ,1985,) ; (g) Utilizin ^ozokl J 19851 ): (h) iChen ot J1 )2005D 
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Table 6. Free-Free Contribution to the 24 yuni Excess 



(1) 



[24] Flux [24] error" 

(mjy) (mJy) 

(2) (3) 



Excess flux at [24] Free-free component 

(mJy) (mJy) 

(4) (5) 



free-free cr excess Excess Type 



(B) 



(7) 



(8) 



GOll. 2691+00. 4208 
G014. 4239-00. 7657 
G027. 0268+00. 7224 
G036. 8722-00. 4112 
G047. 3677+00. 6199 
G047. 4523+00. 5132 
G051. 6491-00. 1182 
G299. 1677-00.3922 
G299. 7090-00. 9704 
G300. 0992-00. 0627 
G305. 4232-00. 8229 
G307. 9784-00. 7148 
G310. 5420+00. 4120 
G311. 0099+00. 4156 
G311. 6185+00. 2469 
G314. 3136-00. 6977 
G321. 7868+00. 4102 

G339. 4392-00. 7791''''' 
0339.7415-00.1904' 

0340.0517+00.6687*' 



18.4 


2.7 


13.1 


71.9 


13.7 


71.2 


17.7 


2.1 


8.7 


21.7 


2.2 


16.3 


44.6 


4.5 


38.9 


8.5 


1.1 


6.3 


17.8 


2.8 


12.7 


7.0 


0.8 


5.4 


5.1 


0.5 


3.1 


7.3 


0.8 


6.1 


16.7 


1.8 


14.2 


21.5 


2.2 


16.1 


5.6 


0.7 


3.3 


149.4 


17.9 


140.9 


8.6 


1.5 


7.4 


22.8 


2.3 


17.8 


12.6 


1.3 


10.7 


7.7 


1.5 




5.4 


0.8 


4.4 


31.6 


3.4 


26.9 



6.0 


46±21% 


2.6 


I 






5.2 


III 


8.5 


98+24% 





I 


7.4 


45+13% 


4.0 


11 


12.1 


31±12% 


6.0 


11 






5.5 


III 


6.1 


48±22% 


2.4 


I 


5.6 


104±15% 





I 


1.9 


61±16% 


2.4 


I 


1.7 


28+13% 


5.5 


II 


6.1 


36+13% 


5.1 


II 


7.0 


43+14% 


4.1 


II 


2.3 


70±21% 


1.4 


I 






7.9 


III 


1.2 


16±20% 


4.1 


II 


4.0 


22±13% 


6.0 


II 


6.0 


47+12% 


4.4 


II 






5.5 


III 


14.4 


54+14% 


3.7 


II 



Note. — (7) Determined by removing the free-free component from the [24] excess and dividing by the [24] error with 
10% calibration uncertainty error added in quadrature with measured error, (a) [24] flux measurement error with 10% 

calibration uncertainty error added in quadrature with measured error, (b) found to have no real excess at [24]. (c) New 
measurement, not in jUzpen et al.||2007| ). 
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Table 7. Stellar and Circumstellar Disk Parameters 



ID 


K-8 


8-24 


Disk Temperature 


^ 




[mag] 


[mag] 


(K) 




(1) 


(2) 


(3) 


(4) 


(5) 


G007.2369+01.4894 






493lt 


„ , . + 0.002 
"■14_o 002 


G008. 3752-03. 6697 






566t^ 


"■2"_0.004 


G009. 4973+05. 2441 










GOll. 2691+00. 4208 


0.34 


1.01 






G014.4239-00.7657 


0.47 


4.78 


191+^ 


o.oi3_o 00, 


G025. 6122+01. 3020 










G027.0268+00.7224 


0.34 


0.56 






G036. 8722-00. 4112 


0.27 


1.26 


427t- 


00072+""'""'^ 
O.U0072_o 00006 


G047.3677+00.6199 


0.93 


1.24 


783t^^ 


0.0028_o 0001 


G047.4523+00.5132 


0.29 


1.25 


355t«, 


0.0012_o 0001 


G051. 6491-00. 1182 


0.19 


1.12 






G229. 4514+01. 0145 






579+* 


058+°""'- 

U.UOB_o 001 


G257.6236+00.2459 






«4t?o 


„ ,0+0.002 
'^■^^-0.002 


G265. 5536-03. 9951-A 










G265. 5536-03. 9951-B 










G269. 5873-05. 8882 










G299. 1677-00. 3922 


0.35 


1.18 






G299. 7090-00. 9704 


0.18 


0.81 






G300. 0992-00. 0627 


0.56 


1.36 


526+- 


n nn7i +0-0002 

0.0021_o 0002 


G305. 4232-00. 8229 


0.79 


1.27 


675li^ 


n nn-?a+00003 
0.0038_o 0002 


G307.9784-00.7148 


0.40 


1.06 


557lf, 


n 00090+"-"'""'^ 

U.0009U_o. 00008 


G310. 5420+00. 4120 


0.25 


0.70 






G311. 0099+00. 4156 


0.37 


2.79 


315t^ 


n nm7+00001 
0.0027_o 0003 


G311. 6185+00. 2469 


0.28 


1.82 


306+g^ 


n nni 9+0-0001 

0.0012_o 0001 


G314.3136-00.6977 


0.22 


1.45 


328+" 


00096+°-°°°°'' 


G321. 7868+00. 4102 


0.77 


1.25 


556t.il 


n nnoi+00002 
0.0021_o 0001 


G339. 4392-00. 7791'' 










G339. 7415-00. 1904 


0.27 


1.93 


31118 


n nnm+OOOOl 
0.0010_o,oooi 


G340. 0517+00. 6687 


1.11 


1.04 


823+^0 


n niK+o.ool 
0.015_o 001 


G347.3777+04.2010 






4871?, 


„ , 3+0.02 
"■-'•^-0.03 



EW(HQ:corT-) "^eff '^ ^^^^ ^ Evolutionary Status 

(A) (K) (kms-1) 

(6) (7) (8) (9) 



-26.15+0.15" 
-0.17+0.08" 
5.11+0.02" 
-13.86+0.16 
4.45+0.03 
6.42+0.09 
-13.89+0.43 
-15.07+0.15 
-22.47+0.24 
6.49+0.02 
-2.89+0.15 
-7.82+0.45" 
9.41+0.06" 
4.93+0.11" 
3.77+0.06" 
2.12+0.03" 
-13.99+0.12 
-10.12+0.16 
-21.65+0.32 
-21.98+0.42 
-13.27+0.69 
-9.97+0.09 
5.21+0.03 
-10.90+0.23 
-12.63+0.29 
-26.81+0.53 
-11.20+0.76 

5.50+0.03 
-39.38+0.08 
-0.02+0.10" 



10350+230 
7260+160 
8320+130 
12350+520 
16900+800 
12520+510 
10510+300 
11740+260 
15700+600 
12260+510 
13100+530 
10220+390 
11530+300 
7660+140 
7040+170 

16450+1030 
12520+505 
13510+510 
13520+510 
13220+520 
12430+510 
9800+130 
10540+255 
9870+130 
12340+505 
8800+150 
12220+510 
14080+520 



145+15 

35±5 
140+15 
220+15 
165+10 
290+20 

25+5 
295+20 
310+25 
150+10 
270+20 
295+20 
240+25 
25+5 
40+5 
10+5 
130+10 
175+15 
255+20 
230+20 
240+20 
270+20 
235+25 
265+20 
205+20 
230+20 

40+5 

240+15 

230+15 

15+5 



H 

TT 



T 

T/B. 

D 



T/Be 
T/Be 
T/Be 

Be 

D 

T 

T 

H 

D 
H 

TT/T 



Note. — (9) H: Herbig AeBc star, TT; T-Tauri star, T; Transitional disk, D: Debris disk. Be: Execss due to free-free emission 
and tliercforc a classical Be-type star, ?: additional measurements are needed to clas sify the nature of th e IR excess, (a) This value 
is EW(Hck), uncorrected for underlying stellar absorption, (b) The [8.0] excess trom lUzpen et al.l 1200711 is not confirmed at [24]. 
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Table 8. Classification Based on lAndrillat et al.l (Il990l ) Spectral Features 



ID Criterion 1 Criterion 2 Criterion 3 Classieal Be 

(1) (2) (3) (4) (5) 



G007.2369+01.4894 
G027.026S+00.7224 
G047.3677+00.6199 
G229. 4514+01. 0145 
G299. 1677-00. 3922 
G300. 0992-00. 0627 
G305. 4232-00. 8229 
G307.9784-00.7148 
G321. 7868-1-00.4102 
G340. 0517-1-00. 6687 



n 


n 


y 


y 


y 


N/A 


N/A 


N/A 


n 


n 


n 


y 


N/A 


y 


N/A 


N/A 


n 


n 


N/A 


n 


n 


N/A 


n 


n 


n 


y 


n 


n 


y 


n 



Note. — (2) Criterion 1: Stars with Ca II in emission and Paschen 
in emission are denoted by y. Stars with Ca II in emission without 
Paschen emission are denoted by n. Stars without Ca II emission 
are denoted by N/A. (3) Criterion 2: Stars with O I in emission and 
Paschen in emission are denoted by y. Stars w^ith O I in emission and 
Paschen in absorption are denoted by n. Stars w^ithout O I in emis- 
sion are denoted by N/A. (4) Criterion 3: Stars w^ith [8.0] excess > 
0.55 magnitudes w^ith Ca II in emission are denoted by y. Stars w^ith 
[8.0] excess > 0.55 magnitudes w^ithout Ca II are denoted by n. Stars 
without [8.0] excess > 0.55 magnitudes are denoted by N/A. (5) Stars 
with criterion 1, 2, and/or 3 with y only are consistent with classical 
Be stars. 
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Fig. 1. — (u pper left) The relati on between projected rotational velocity and He I A 4471 
A FWHM of iMunari et al.l ( 20051 ) stellar models, where the multiple points at each velocity 
reflect measurements over a range of effective temperatures. The He I A 4471 A line contains 
Stark broadening that results in a greater dispersion in the FWHM measurement. The 
solid line shows the linear least squares fit. (upper right) The relation between projected 
rotational velocity and Mg II A 4481 A FWHM. (lower left) The relation between projected 
rotational velocity an d Fe I A 44 7 6 A F WH M. (lower right) A com parison of projected 
rotational velocities of iRoyer et al.l (l2002al ) and ISlettebak et al.l (119751 ) versus our rotational 
velocities. The solid li ne shows the 1:1 rel ationship in this panel. Our m easurements are 
consis tent with those of iRoyer et al.l (l2002al ) but ~10% larger than those of ISlettebak et al. 
( 119751 ). as explained in the text. 
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Fig. 2. — (upper left) The relation between model ratios of EW Mg II A 4481 A/EW He I 
A 4009 A and effective temperature. The multiple data points at each abscissa/ordinate are 
due to multiple measurements of the same model and shows the systematic errors associated 
with this relationship. The solid curve shows the best fit relation of EW ratio to effective 
temperature for the models with 150 km s~^ projected rotational velocity, (upper right) The 
relation between EW Mg II A 4481 A/EW He I A 4471 A and effective temperature, (lower 
left) The relation between EW Mg II A 4481 A/EW Fe I A 4233 A and effective temperature. 
This relationship depends upon rotational velocity, and the dotted curve shows the best 
fit relation of EW ratio to effective temperature at 250 km s"^. (lower right) The relation 
between EW Mg II A 4481 A/EW Ca I A 4227 A and effective temperature. 
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Fig. 3. — He I A 5876 A profiles (d enoted by the dashe d vertical line) of mid-infrared ex- 
cess stars overplotted with best fit iMunari et al.l (120051 ) stellar models. Red curves show 
the stellar models for stars in which the measured stellar temperature falls in between the 
two nearest stellar models. Blue curves show the stellar models for stars in which the mea- 
sured temperature does not fall in between the nearest stellar models. G014. 4239-00. 7657, 
G229.4514+01.0145, G265.5536-03.9951-A and B are overplotted by blue curves indicating 
the discrepancy between the He I line and measured temperature. 
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Fig. 4. — H e I 58 76 A profiles of mid- infrared excess stars overplotted with best fit 
Munari et al.l (120051 ) stellar models using the same symbolism as Figure [31 G299.1677- 
00.3922 and G307. 9784-00. 7148 are overplotted by blue curves indicating the discrepancy 
between the He I line and measured temperature. 
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Fig. 5. — Ha profiles of mid- infrared excess stars. Tlie majority of our sources exliibit a 
double-peaked emission profile. 
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Fig. 6. — Ha profiles of mid- infrared excess stars. Tlie majority of our sources exliibit a 
double-peaked emission profile. 
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Fig. 7. — The spectral region 8400-8600 A showing several Paschen lines, Pal5-19, Ca II A 
8498 and 8542 A and O I A 8446 A. 
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Fig. 8. — The spectral region 8400-8600 A showing several Paschen lines, Pal5-19, Ca II A 
8498 and 8542 A and O I A 8446 A. 
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Fig. 9. — The spectral region 8400-8700 A showing the three stars within our sample that 
exhibit Ca II emission. iMunari et al.l (120051) stellar models of nearest measured temperature 
are shown below. All three sources also exhibit O I emission, and G027.0268+00.7224 
exhibits Paschen emission. 
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Fig. 10. — Fe II A 5317 A profiles of mid-infrared excess stars. Fe II is indicated by the 
dashed vertical line. Fe II emission is not very common among our sources. 
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Fig. 11. — Fe II A 5317 A profiles of mid-infrared excess stars. Fe II is indicated by the 
dashed vertical line. Fe II emission is not very common among our sources. 
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Fig. 12. — HR diagram for sixteen of our target stars overlaid with lSiess et al.l (120001 ) stellar 
evolutionary models for 7,000 < Tg// < 17,000 K. The symbols show the various methods 
used to determine the absolute magnitude, diamonds for Hipparcos parallax, triangle for 
Tycho parallax, and asterisks for stars with absolute magnitudes taken from the literature. 
The zero-age main-sequence is shown by the solid line. Isochrones are shown by the dotted, 
dashed, and dash-dot relations. The majority of our stars lie above the main-sequence and 
are consistent with 3-5 Mq stars. The three coolest stars are above the main-sequence and 
consistent with pre-main-sequence objects. 
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Fig. 1 3. — A comparison of pr ojected rotationa l veloc ities of B8-9.5 V stars from lAbt et al.l 
(I2OO2I ). classical Be stars from lHanuschik et al.l (119961 ) . and our GLIMPSE sample. U t ilizing 
a KS test we find that there is no signi fi cant p robability, < 1% for the lAbt et al.l (|2002| ) 
sample and 8% for the iHanuschik et al.l (Il996l ) sample, that this sample may be derived 
from either population. 
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Fig. 14. — A i^-8 /iin vs. 8-24 /xni color-color plot of our objects {open symbols) and classical 
Be stars (asterisks) fromlJaschek fc Egretl (119821 ) with MSX A band and IRAS 25 //ni (from 
Zhang et al.ll2004j . |2005| ) were utilized as 8 /xm and 24 ^m measurements respectively. The 
arrow shows the reddening vector for Ay=5.0. The solid curve denotes a B8V main-sequence 
with Ay =1.0 and increasing contribution of free-free flux at 8 yuni. The solid symbols denote 
5%, 33%, and 50% circumstellar contributions to the total flux at 8 /im. The dashed curves 
denote B8V main-sequence stars with Ay =1.0 and increasing contributions from 800 K 
and 300 K blackbody circumstellar components, respectively. This plot shows that excess 
emission owing to free-free processes can be distinguished from cool dust (< 300 K). Systems 
containing the former lie to the left of the solid line while the systems with the latter lie to 
the right. In this color-color space, free-free and hot dust excesses are indistinguishable. 
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Fig. 15. — Spectral energy distribution of one of our target stars, G036. 8722-00. 4112. The 
thick sohd hne is the distance-normahzed stellar model. The thin solid line is the stellar 



model with extinction and circumstellar blackbody component from lUzpen et al.l (120071 ). 
The dash-dot curve is the modeled free-free component normalized using the Ha flux. The 
dashed-dot-dot-dot curve is the sum of the stellar model and the free-free component. The 
free-free component is not sufficient to explain the measured excess at 24 fim. The residuals 
are for the stellar model with circumstellar blackbody component. 
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Fig. 16. — (a) EW(HQ;corr) VS. disk temperature. There is a strong anti-correlation between 
these two parameters, (b) Ha equivalent width vs. stellar effective temperature. There is 
no correlation between these two parameters for Type I, Type II, or Type III sources. 
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Fig. 17. — An inter-comparison of -j^, age, and Ha equivalent width. Using the younger 
age, when available, age is correlated with Ha equivalent width with a correlation coefficient 
of 0.97 indicating > 99% probability that the two parameters are correlated for these seven 
stars. Using the older age when available, there is no correlation between age and -^ or 
between Ha equivalent width and -j^. 
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Fig. 18. — A flow chart showing how we determined the type of circumstellar disk system 
using echelle spectra and an estimate of the contribution from free-free emission for 19 
GLIMPSE stars. Two MSX stars are characterized: G007.2369+01.4894, where extensive 
data exists in the hterature, and G229. 4514+01. 0145, which exhibits spectral features similar 
to G007.2369+01.4894. 



